Microscopic dynamics in liquid metals: the experimental point of view. 



Tullio Scopigno^ and Giancarlo Ruocco 1 

1 INFM CRS-SOFT and Dipartimento di Fisica, Universita di Roma " LaSapienza" , 1-00185, Roma, Italy 
Francesco Sette 2 

2 European Synchrotron Radiation Facility, B.P. 220 F-38043 Grenoble, Cedex France. 
(Dated: February 2, 2008) 

The experimental results relevant for the understanding of the microscopic dynamics in liquid 
metals are reviewed, with special regards to the ones achieved in the last two decades. Inelastic 
Neutron Scattering played a major role since the development of neutron facilities in the six- 
ties. The last ten years, however, saw the development of third generation radiation sources, 
which opened the possibility of performing Inelastic Scattering with X rays, thus disclosing previ- 
ously unaccessible energy-momentum regions. The purely coherent response of X rays, moreover, 
combined with the mixed coherent/incoherent response typical of neutron scattering, provides 
enormous potentialities to disentangle aspects related to the collectivity of motion from the single 
particle dynamics. 

If the last twenty years saw major experimental developments, on the theoretical side fresh ideas 
came up to the side of the most traditional and established theories. Beside the raw experimental 
results, therefore, we review models and theoretical approaches for the description of microscopic 
dynamics over different length-scales, from the hydrodynamic region down to the single particle 
regime, walking the perilous and sometimes uncharted path of the generalized hydrodynamics 
extension. Approaches peculiar of conductive systems, based on the ionic plasma theory, are also 
considered, as well as kinetic and mode coupling theory applied to hard sphere systems, which turn 
out to mimic with remarkable detail the atomic dynamics of liquid metals. Finally, cutting edges 
issues and open problems, such as the ultimate origin of the anomalous acoustic dispersion or the 
relevance of transport properties of a conductive systems in ruling the ionic dynamic structure 
factor are discussed. 
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I. INTRODUCTION 



Liquid metals are an outstanding example of systems 
combining great relevance in both industrial applica- 
tions and basic science. On the one hand they find 
broad technological application ranging from the produc- 
tion of industrial coatings (walls of refinery coker, drill 
pipe for oil search) to medical equipments (reconstruc- 
tive devices, surgical blades) or high performance sport- 
ing goods. Most metallic materials, indeed, need to be 
refined in the molten state before being manufactured. 

On the other hand liquid metals, in particular the 
monoatomic ones, have been recognized since long to be 
the prototype of simple liquids, in the sense that they 
encompass most of the physical properties of real fluids 
without the co mplications which may be pr esent in a par- 
ticular system l)Balucani and ZoppHll983(l . 

In addition to that, metallic fluids such as molten 
sodium, having similar density and viscosity as water, 
find application as coolant in nuclear reactors. 

The thermodynamic description of liquid metals can 
be simplified by assuming a few parameters. Usually, if 
compound formation is weak physical theory alone can 
be used while, if there is strong compound formation, 
chemical theory alone is used. The lowest-melting liquid 
metals are those that contain heavier elements, and this 
may be due to an increase in ease of creating a free- 
electron solution. Alkali metals are characterized by low 
melting points, and they tend to follow trends. Binary 
associating liquids show a sharp melting point, with the 
most noticeable example being mercury (T m = 234 K). 
Melting points can be lowered by introducing impurities 
into the metal. Often, to this purpose, another metal 
with a low melting point is used. Mixing different metals 
may often result in a solution that is eutectic. In other 
words, from Henry's law it is understood that a melting 
point depression occurs, and the system becomes more 
disordered as a result of the perturbation to the lattice. 
This is the case, for instance, of the well known eutectic 
Pb-Tin alloy, widely used in soldering applications (T m = 
453 K). 

Until the sixties the understanding of the physical 
properties of metals proceeded rather slowly. It was 
John Ziman, indeed, who made t he theory of li quid met- 
als respectable for the first time l)Zimanlll960(l . and the 



Faber-Ziman theory, developed in 1961-63 and dealing 
with electronic and transport properties, is attractively 
introduced in Faber's book, which is an ex cellent treatis e 
of the physical properties of liquid metals (jFaberl fl972^l . 

The other t ext wh i ch ca n be considered a classic is 
March's book llMarcbi ll968h along with the more recent 
ijMarcbl Il99d|) , which provides a comprehensive overview 
over liquid metals. It is from these texts that a first clear 
definition of liquid metal can be outlined. At first glance, 
indeed, the words " liquid metal" are self-explanatory: by 
definition any metal heated to its melting point can be 
cast in this category. Liquid metals, however, are implic- 
itly understood to be less general than the above defini- 
tion, and no literature clearly states an exact definition. 
Although no precise agreement has been made, there are 
certain characteristics shared by liquid metals, descend- 
ing from a close interplay between ionic structure, elec- 
tronic states and transport properties. 



The book of Shimoji l|ShimoiiL Il977|) deals with the 
fundamentals of liquid metals in an elementary way, cov- 
ering the developments achieved after the first book by 
March. It does not address, however, the dynamical 
properties in great detail. 



Addison's book iAddisonl lillStl is much like March's 
general book, but is more focused on applications of al- 
kali metals, especially on their use in organic chemistry. 
In addition, Addison discusses many methods for puri- 
fying and working with liquid alkali metals. March is 
more theoretical whereas Addison is practical, but both 
authors focus on a thermodynamic explanation of liquid 
metals. 

For an appealing general introduction to the physics 
and chemistry of the liquid-vapor phase transition (be- 
yond the scope of t his review) the reader shou ld certainly 
make reference to l|Hensel and Warren! 1 1999(1 . which also 
provides a bird eye view of the practical applications of 
fluid metals, such as high-temperature working fluid or 
key ingredients for semiconductor manufacturing. 

There are, then, a number of books which 
are more general and more specific at the same 
time, in the sense that they deal with with the 
wider class of simple liquids (including noble flu- 
ids, hard sphere fluids etc.), but they are mainly 
concerned with structural and dynamical proper- 



ties only llBalucani and ZoppiL [1983; 


Boon and Yipl 


198(i lEs-elstafl 


119671 




Hansen and McDonald 


eh 


1986: 


March and Tosil 1 


19911) 


They are practically in 


idible 



for those aiming at a rigorous approach to the statistical 
mechanics description of the liquid state. 

It can be difficult to find an exhaustive updated 
database of the physical properties of liquid metals, es- 
pecially a s far as dynamics is con cerne d. But the hand- 
books of l|Iida and Guthriel Il993[) and l|Ohsd. Il985[) are 
remarkable exceptions, with the second one specifically 
addressing liquid alkali metals. 
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A. Historical background 

Early phenomenological approaches to the study of re- 
laxation dynamic in flui ds can be dated back to the e nd 
of the nineteen century l|Kelvinl Il875l iMaxwell Il867|) . 

Only in the mid twentieth century, however, it was re- 
alized that a deeper understanding of the physical prop- 
erties of liquids could have been reached only through 
a microscopic description of the atomic dynamics. This 
became possible through the achievements of statistical 
mechanics which provided the necessary tools, such as 
correlation functions, integral equations etc. The math- 
ematical difficulties related to the treatment of real liq- 
uids brought to the general attention the importance of 
simple liquids, as systems endowed with the rich basic 
phenomenology of liquids but without the complications 
arising, for instance, by orientational and vibrational de- 
grees of freedom. 

As a consequence, the end of the fifties saw major ex- 
perimental efforts related to the development of Inelastic 
Neutron Scattering (INS) facilities which, as we shall see, 
constitutes a privileged probe to access the microscopic 
dynamics in condensed mat ter and, in particular, in the 
liquid state ijEgelstafll H959) . A sizable library of exper- 
imental data on liquid metals has been constituted since 
then, realizing the prototypical structural and dynamical 
properties of these systems, representative of the whole 
class of liquids. 

In the sixties, the advent and the broad diffusion of 
computational facilities brought a new era for two main 
reasons: on the one side, rea listic compu ter simulation 
experiments become possible l|Schift Il969|) . on the other 
side the new computation capabilities greatly facilitated 
the interpretation of INS experiments. For instance, new 
protocols for accurate estimates of the multiple scattering 
contribution a ffecting neutron scattering were proposed 
l)Coplevl Il975|) . The theoretical framework of the Inelas- 
tic Neutron Scattering, and the guidelines to interpret 
the results, have been rev i ewed in the classical textbooks 
llLovesevi Il987l iMarshalli Il97lfl . 

The dynamics of liquid metals has been extensively 
investigated by INS and computer simulations with the 
main purpose of ascertaining the role of the mechanisms 
underlying both collective and single-particle motions at 
the microscopic level. In the special case of collective 
density fluctuations, after the semin al inelastic neutron 
scatt ering study by Copley and Rowe ijCoplev and Rowel 
Il974|) an d the famous mo lecular dynamics simulation of 
Rahman ijRahmanl 11974^ in liquid rubidium, the inter- 
est in performing more and more accurate experiments is 
continuously renewed: it was soon realized, indeed, that 
well-defined oscillatory modes could be supported even 
outside the strict hydrodynamic region. In molten alkali 
metals, moreover, this feature is found to persist down to 
wavelengths of one or two interparticle distances, mak- 
ing these systems excellent candidates to test the various 
theoretical approaches developed so far for the micrody- 
namics of the liquid state. 



Up to ten years ago the only experimental probe ap- 
propriate to access the atomic dynamics over the inter- 
particle distance region were thermal neutrons, and using 
this probe fundamental results have been gained. There 
are, however, certain limitations of this technique which 
can restrict its applicability: First, the presence of an 
incoherent contribution to the total neutron scattering 
cross section. If on one hand this allow to gather a richer 
information, being simultaneously sensitive to collective 
and single particle dynamics, on the other hand poses 
the problem of decoupling the two contributions, when 
aiming at the determination of collective properties only 
(i.e. of the coherent dynamic structure factor S(Q,uj)). 
In liquid sodium, for instance, the incoherent cross sec- 
tion dominates; even in more favorable cases (Li, K) at 
small Q the intensity of the collective contribution is low, 
and its extraction requires a detailed knowledge of the 
single particle dynamics. 

The second reason is dictated by the need of satisfy- 
ing both the energy and momentum conservation laws 
which define the (Q — E) r egion accessible to the probe 
l|Balucani and ZoppiL ^)83). Roughly speaking, when 
the sound speed of the system exceeds the velocity of 
the probing neutrons (~ 1500 m/s for thermal neutrons) 
collective excitations can hardly be detected for Q values 
below, let's say, Q m , the position of the main diffraction 
peak of the sample, which is the region where collec- 
tive properties show the richer phenomenology. As we 
shall see in section Ull. A. 51 given a certain kinematic re- 
gion accessible to neutrons (basically ruled by their ther- 
mal energy), by virtue of the m -1 / 2 dependence of the 
sound velocity of an atomic system, the higher its atomic 
number, the wider is the accessible energy-momentum 
region of the excitations which can be studied. Tak- 
ing as an example alkali metals, indeed, accurate INS 
data are availabl e for heavier e l ements such as rubidium 
(v ~ 1260 m/s) llChieux et all Il996t ICoplev and Rowel 
11974 iPasaualini et all Il999l) and cesium (v ~ 970 m/s) 
l|Bodensteiner etalL Il992|) . while more difficulties are 
met in the case of lighter atoms. In particular, lithium 
represents the most critical case due to its high sound 
speed (v ~ 4500 m/s) and to the weak scattering cross 
section which, moreover, results from comparable values 
of the coherent and incoherent contributions: for this 
reason INS aiming to the study of colle ctive properties 
of Li represented a v ery hard challenge l)De Jond . Il993l 
iDe Jong et afl.ll993fl . From a general point of view the 
main outcome of most of these early INS experiments, as 
far as collective properties are concerned, is the evidence 
of inelastic excitations in S(Q, to) which have been neces- 
sarily analyzed within simpl e models such as the da mped 
harmonic oscillator (DHO) l)Fak and Dornei(ll997|) . suit- 
able to extract reliable and resolution-corrected infor- 
mation on the peak positions but not about the de- 
tail of the whole lineshape. Some additional informa- 
tion ha ve been achieved, for inst ance, in the case of 
cesium l|Bodensteiner et all Il992|) . where information 
about an average relaxation time have been extracted 
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utilizing Lovesey's viscoelastic model (jLovesevMl97l|) or, 
more recently, in molten potassium, where a general- 
ized hydrodynamic treatment as the one d escribed in sec. 
III.FI is undertaken iCabrillo et all l2002ft and electronic 
screening effects h ave been explicitly taken into account 
llBove <l nl\.Uml\. 

Paralleling the development of INS facilities, new 
ideas arose on both theoretical and numerical fields 
from 1975 and in the intervening deca des, driven by ki- 
netic theory appli e d to Enskog's fluid (|Ane_Y_^iid Alder , 
'I983t lAllev et all Il983t Ide Schepper and Cohenl I198C : 



tering iBurke]Lll99l became real thanks to the advent of 
the third generation sources ijMasciovecchio et"q/lll996|) . 



Ide Schepper et aU. Il983t iKamgar-Parsi et ffla il987t> . al 

lowing to describe the hydrodynamic region in terms of 
three pole approximation, or to reproduce the dynamic 
structure factor at wavelengths comparable to the in- 
verse mean particle distan ce in terms of extended heat 
mode l|Cohen et all Il987h . Kinetic approaches where 



and by Mode Coupling theorv llBengtzelius et all 


1984 




1976 


Sioarcn, 1979, 1980: Sioarcn and Siolandcr, 1978, 


1979 


Wahnstrom and SioarenL |l982) . 



Turning the attention to numerical advances in the 
liquid metals field, the major achievement are proba- 
bly relat ed to the introduction of the pseudo-po t ential 
concept llAshcroft and StroudL fl978t I Austin et qZl Il962t 
iHarrisonl Il966t iHeind . Il970(l which, beside offering a 
deeper comprehension of physical properties such as elec- 
trical resistivity, provided a clue for realistic numerical 
simulations. In Molecular Dynamics, indeed, the choice 
of a realistic interatomic potential - i.e. a potential model 
able to reproduce structural properties - is crucial to al- 
low the determination of the dynamics of the system via 
the integration of the classical Newton equations. Ex- 
ploiting the pseudo-potential theory, it has been possible 
to express the atomic interaction as a sum of pairwise in- 
teractions, ruled by an effective density dependent inter- 
action. In this respect, one of the most successful expres- 
sions is the effectiv e potential proposed for alkali metals 
l|Price et all Il970|) . 

The numerical simulation framework is particularly 
useful since the single particle and the collective dynam- 
ics can easily be investigated within technical restrictions 
due to the finite box size (defining the minimum acces- 
sible wavevector) and computation time (related to the 
statistical quality and to the energy resolution of the cal- 
culated spectra). Broadly speaking the features of the 
atomic collective motion i.e. the details of S(Q, ui) line- 
shape, as outcome of MD run, turns out to be less noisy 
and more straightforward than the correspondent INS 
results: no absolute normalization is required, no mix- 
ing between coherent/incoherent dynamics occurs and, 
above all, basically no resolution corrections are needed. 

The major experimental breakthrough, however, hap- 
pened in the last ten years when x-rays came up by the 
side of neutrons to study the collective dynamics in a 
similar frequency and wavelength region. The intriguing 
theoretical possibility of performing Inelastic X ray scat- 



disordered systems llRuocco et all 


1996t IScopigno et all 


2003: ISette et all Il998l 


ISinn et al. 


, 2003). In this case, 



the cross section is mainly coherent, and the combination 
of the two techniques can in principle serve to disentangle 
the two contributions. Unfortunately, such complemen- 
tarity has not yet been exploited at full. 



B. Dynamical aspects 

There are mainly two routes to approach the dynamics 
of a viscous melt. The first one stems from a quasi- 
crystalline picture, and relies on the observation that, 
often, the diffusion coefficient D linearly depends on 1/T. 
The same dependence, in fact, is induced in crystalline 
solids by vacancies and defects. This analogy suggests 
that diffusion in liquids is an activated process, and many 
attempts have been done to relate the activation energy 
in the liquid to the thermodynamics of its solid. 

The other point of view is the kinetic theory, a gas-like 
picture where the correlation functions are different in 
view of the density which is typically much higher than 
in the gas state. Within this framework an expression 
for the behavior of diffusion coefficient and viscosity can 
be gained in terms of the friction coefficient £. When a 
particle of the melt is moving with constant velocity v a 
net retarding force results from the different rate between 
front and back collisions of the form F = —£,v. For an 
hard sphere gas it turns out from Fick's law of diffusion 
that £ = 8 " g g tT -*' T y/nmkBT, where p is the atomic den- 
sity, m the atom mass, and 4ira 2 g(<r) the density prob- 
ability of finding two units at distance a. The diffu- 
sion coefficient is therefore D — while the viscosity 

1) = gf^: l)Longuet-Higgins and Popjel Il956h . These sim- 
ple expressions turn out to describe remarkably well the 
dynamics of liquid metals as long as the one of other 
simple fluids. However, better quantitative agreement 
can be obtained introducing the velocity autocorrelation 
function ^p(At) — ^ t+ ^2) , whose time integral de- 
termines the diffusion coefficient. The i/j(t) plays a cen- 
tral role in liquid dynamics, not only for providing a rig- 
orous way to calculate D, but also because through its 
Fourier transform one can grasp an insight into the detail 
of the interatomic interactions. It can be accessed either 
by molecular simulations (the calculation trivially follows 
from its definition) or experimentally, mainly by Inelastic 
Neutron Scattering (INS) with the methods detailed in 
section Mil 

Broadly speaking, ip(At) is related to the knowledge 
of G s (r, t) i.e. the probability that a given particle trav- 
els a distance r in the time interval t. INS is always 
sensitive to a combination of G s (r, t) and Gd(r, t), being 
this latter quantity the probability of finding two distinct 
particles at a space distance r and time distance t. The 
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way to separate these two contributions mixed in the in- 
strumental response is one of the major conundrum of 
the neutron scattering technique. The direct knowledge 
of the coherent response, gained by means of IXS, opens 
the possibility to the understanding of the high frequency 
modes that have been seen to survive since the famous 
INS experi ments on the highly coh erent scatterer liquid 
Rubidium (C opley and Rowel Il974j) . In particular, from 
a solid point of view, the main issue is to ascertain the 
relation between these modes and the phonon excitations 
in the corresponding polycrystal just below the melting. 
In a perfect crystal (i.e. a periodic assembly of atoms 
or ions), indeed, atomic dynamics is mainly vibrational, 
characterized by normal modes which are plane waves, 
due to the harmonic nature of the interatomic forces. 
Consequently, a well defined ratio (the sound velocity) 
exists between the frequency and the wavevector of the 
density fluctuations, named in this case phonon. In real 
crystals, therefore, the energy spectrum of a scattered 
probe results in sharp peaks, whose lincwidth is related 
to the presence of anharmonicities or lattice defects. In a 
liquid things are more involved: other effects arise besides 
anharmonicity, both structural (the average atomic posi- 
tions are randomly distributed) and dynamical (mass dif- 
fusion and activated processes join the purely vibrational 
motion) . But even in this complex scenario distinct peaks 
survive and one can extract dispersion relations. 

Interestingly, the relation between crystals and liquids 
involves more than the mere similarity between the sound 
velocity of the crystalline acoustic branches and of the 
mode observed in the liquid (i.e. the low Q limit of 
the dispersion relation). Despite the lack of periodic- 
ity intrinsic to the inherent liquid structure, indeed, the 
presence of some residual correlation (testified by the os- 
cillation in static structure factor) seems to warrant a 
support for the existence of umklapp modes similar to 
the one existing in crystals. Such processes, i.e. the 
presence of inelastic modes characterized by wavevectors 
which differ by multiples of the reciprocal lattice spac- 
ing, have been early reported by means of INS in liquid 
lead llCocking and Egelstafli Il965at iDorner et oil Il965t 
lRandolpTHmcH5mlrwT^)6(fi and more recently in liquid 
lithium by means of IXS iScoDigno et all l2001bl> . As 
correctly pointed out by Faber, the presence of these ex- 
citations does not imply the existence of genuine high Q 
modes, it rather indicates that um klapp proces ses may 
occur in liquid as much as in solid ilFabeitll972ft . 

From the liquid point of view the interest in this phe- 
nomenology lies in the challenging extension of the sim- 
ple hydrodynamics, describing the density fluctuations in 
the long wavelenght limit, down to the lcnghtscale of the 
mean interparticle distances. As it will be shown in the 
following, such an extension relies on serious and some- 
times not fully justified assumptions necessary to walk 
in the uncharted and perilous territory between hydro- 
dynamics and single particle regimes. 



C. Peculiarities of liquid metals 

Apparently at odds with the previously mentioned 
classification of liquid metals as prototype of simple liq- 
uids, even in the simplest monoatomic case, metallic flu- 
ids are actually two component systems. The interplay 
between electron and ions, indeed, is an intrinsic aspect of 
liquid metals, and a rigorous approach should therefore 
mimic the formalism utilized for binary mixtures. For 
many aspect, however, one might be interested in ionic 
properties only, and as far as atomic dynamics is con- 
cerned, this seems to be the case. In such circumstance, 
one can look at a liquid metal as an ionic assembly whose 
interaction is mediated by the conduction electron gas. 
The treatment is in such way reduced to a one component 
system, like for noble fluids one can introduce a pairwise 
interaction, but this latter will be ruled by a density de- 
pendent pseudopotential. 

Although within the pseudopotential approach many 
results for liquid metals are qualitatively similar to those 
for ordinary non conductive fluids, some remarkable dif- 
ferences exist. One of the most relevant of these dif- 
ferences concerns peculiar structural properties involv- 
ing short range order: in several liquid metals the static 
structure factor exhibits an asymmetry or even a shoul- 
der just above the main peak. The origin of this anomaly 
has been highly debated, and ascribed to the peculiar 
shape of the interaction potential in those metals in whic h 
the hard sphere description fails l)Tsav and Wangl 11994^1 . 
More specifically, it can be interpreted in terms of a re- 
pulsive interaction composed by an hard sphere part plus 
an adjacent ledge induced by electronic effects, by a cur- 
vature change occurring at the nearest neighbor distance, 
or by the interplay with Friedl oscillations. 

The reported peculiarities extend also to the dy- 
namics: while in the long wavelength limit they 
are expected to behave similarly to non conduc- 
tive fluids, at finite wavevectors their departure 
from ordinary hydrodynamics can be in principle in- 
fluenced b y the high values of the therma l con- 
ductivity llFaberl Il973 IScopigno and Ruoccol 12004 
ISingh and Tankeshwait 120031 12004^ One of the most 
striking quantitative differences with ordinary fluids con- 
cerns the "visibility" of the inelastic features, i.e. the 
inelastic to elastic ratio, which seems t o be related to 
the softness of the interaction p otential ijBalucani et all 
1992[ lCanales and Padrolll999j) . This latter is responsi- 
ble, for example, for the very favorable inelastic to elastic 
ratio which makes alkali metals ideal systems to study 
collective properties. 

Though the whole dynamics of liquid metals seems 
to be conveniently rationalized treating them as ordi- 
nary fluids interacting via an effective, density depen- 
dent, pairwise interaction potential, there is an alter- 
native route which explicitly takes into account elec- 
tronic screening effect on the ionic dynamics, relying 
on the introduction of a suitable model for the wave- 
length dependent dielectric function. In this way one 
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his able to test different approximations comparing the 
predictions for the mechanical compressibility (or, equiv- 
alen tlv, for the sound velocity) with the e xperimental val- 
ues l)Bove et al I 12003 ISaid et oiJ . l20f)l . 



D. Why this review 

As previously pointed out, some books offer a broad 
coverage of the physics and chemistry of liquid metals, 
but none of them is focused on the dynamical aspects. On 
the other side, books dealing with liquid dynamics do not 
spot on liquid m etals in great detail. Egelstaff's review 
(Egelstaff, 1966) offers an exhaustive coverage of trans- 
port phenomena in liquids but dynamical properties of 
liquid metals ar e marginally addressed. Th e Copley and 
Lovesey review ijCoplev and Lovesevl Il975j) provides an 
invaluable sight on the dynamics of simple liquids, and 
it presents many results obtained for liquid metals, but 
again does not emphasize the peculiarity of these sys- 
tems. Moreover, it also addresses in detail numerical 
results and last, but not least, it is necessarily not up 
to date given the large amount of experiments recently 
performed. The last decade, indeed, has seen important 
advances in the dynamics of liquid metals, especially on 
the experimental side, driven by the advent of new X- 
ray facilities and by the upgrades of the neutron ones. 
On the theoretical side, in the eighthes, impressive re- 
sults have been achieved working on hard sphere models 
and/or mode coupling approaches which, again, remain 
uncovered in the Copley and Lovesey review. To con- 
clude, it is certainly worth m entioning the m ore recent 
review article by P. Verkerk ijVerkerkl 12001). which of- 
fers a clear overview of the theoretical models developed 
sofar for liquid dynamics, presenting a selection of exper- 
imental results for liquified rare gases, molten metals and 
binary mixtures. 

Given this background, it seemed to us helpful to focus 
on the experiments on liquid metals, and to discuss and 
summarize the results and their interpretations in terms 
of the existing theories, trying to emphasizes advantages 
and weakness of each approach. This turned out to be 
a difficult task, given the broadness of the matter, and 
we necessarily had to make some choices. We left out, 
for instance, mixtures and alloys, and we tried to focus 
on the most recent experimental achievements, say of the 
last ten years: in most cases we quickly reference to older 
results, unless they are particularly relevant in view of the 
most recent ones. 

The review is organized as follows. In section [H] we 
present different theoretical approaches to the dynamics 
of liquid metals. We develop in parallel subsections the 
treatment of the self and collective properties which are, 
in turn, organized according to the different wavelength 
domains: hydrodynamic, non-hydrodynamic, single par- 
ticle. We also include a subsection dealing with hard 
sphere treatment and one presenting the ionic plasma 
approach, which is peculiar of conductive systems. In 



section Hill we describe the experimental approach to the 
investigation of microscopic dynamics in liquids, outlin- 
ing the basics of the inelastic scattering problem. Since 
the case of X rays is relatively newer, we decided to treat 
it in detail, but continuous reference is made to neutron 
scattering in an effort to emphasize merit, drawbacks and 
complementarities of the two methods. Section Hvl is the 
bulk part of this paper, in which the experimental results 
are reviewed and ordered element by element. Here, we 
make constant reference to section [n] to recall the differ- 
ent approaches utilized by different authors to describe 
the experimental result. In section finally, we try to 
summarize the arising scenario, pointing out the issues 
which, in our opinion, deserves further investigations and 
trying to draw, when possible, some conclusive pictures. 

II. THEORETICAL BACKGROUND 
A. General overview 

1. Some basic definitions 

The investigation of microscopic dynamics of an en- 
semble of N identical atomic or molecular units usually 
proceeds through the study of correlation functions of 
dynamical variables, i.e. of functions of the phase space 
variables, defined as the 6iV positions r^(i) and momenta 
Pi(t) = mviit) of the particles. Relevant dynamical vari- 
ables are those stemming from the microscopic density 
p(r,t) (whose average is related to the number density 
p = (p(r,t))), momentum density J(r,t) and kinetic en- 
ergy density E(r,t): 

p(r,t) = ^=Yl 6 ( r - r ^ 

v i 

J(r,t) = -^^v i (t)5(r-r i (t)) 
Eir^^-^^mvUmr-r^t)) (1) 

i 

In many cases, the study of the dynamics of a tagged 
particle i can be of interest, and it relies on similar defi- 
nitions for the single particle dynamical variables: 

p s {r,t) =*(r-r,(f)) 
J.(r,t) =Vi(t)«(r -!■<(*)) 

£ s (r,i) = im Wl 2 (^r-r^)) 

The well known van Hove distribution functions 
Gd(r,t) and G s (r,t) are related in the classic (not quan- 
tum) case to the microscopic self and collective densities 
through: 

G(r,t) = G s (r,t)+G d (v,t) 
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with 



Consequently 



G *( r '*) = ^(E d(r + ri(0) ~ r ^* )) ) 



\ i jjti I 

As we shall see, experiments usually give information 
on the correlation functions in the reciprocal, Q, space. 
Therefore, it can be useful to define the space Fourier 
transform of the microscopic quantities previously intro- 
duced: 



p(Q,t) = 

J(Q,i) = 
E(Q,t) = 



—y 



-iQ-ri(t) 



iQTi(t) 



and similarly for the single particle variables. 

The time evolution of these microscopic variables is 
cross-linked and can be studied by replacing them with 
their average over small but statistically significant vol- 
umes. A closed set of equations, basically conservation 
laws and constitutive relations, can easily be written for 
these averages 1 . Once the expression for a(Q, t) is known 
(a = p,3,E) one can calculate 



(2) 



In which a Q (Q) are the appropriate dynamical vari- 
ables, and C is the Liouville operator ruling the time evo- 
lution in the configurational space. The (...)n indicates 
thermal averages evaluated over the iV-particle ensemble 
(from now on we will omit the subscript N). 

In particular, the autocorrelation functions of the mi- 
croscopic density (both self and collective) play a privi- 
leged role, and we define, therefore: 



F(Q,t) = (p(Q,t)p(-Q,0)) 
S(Q) = (\p(Q,0)\ 2 ) 

F s (Q,t) = ( Ps (Q,t) Ps (-Q,0)) 
$ s (Q,t) = F s (Q,t) 



(3) 



1 From here on we will implicitly assume that we are not dealing 
with microscopic quantities but rather with the hydrodynamic 
quantities resulting from their averages 



HQ,t) 



1 V / p iQ-n(°)p-iQ-*iV) 



NS(Q) 



1 J- / e «J-ri(«V* e -<Q-**C°)) (4) 



NS(Q) 
N ^ 



e iQTi(0) e -iQT 4 (t) 



(5) 



The above autocorrelation functions F(Q, t) and F S (Q, t) 
are, in fact, connected through their time- Fourier trans- 
form to the self and collective dynamic structure factor 
S S (Q, oj) and S{Q, uj), respectively. These latter, in turn, 
are the experimentally accessible quantities in neutrons 
and X-rays inelastic scattering experiments and there- 
fore, in the following, we will mainly refer to the density- 
density correlation functions. 



2. Spectral moments 

Before illustrating some models for the evolution of the 
density autocorrelation function, it is worth to recall here 
some basic relations involving the frequency moments of 
the dynamic structure factor. These can be very useful 
to the experimentalist, as a way to normalize the data 
(an example is given in section IIII.BI as well as to any 
theoretical approach, as a direct test of sum rules - for 
example, we shall see how the hydrodynamic expression 
for the density-density time correlation function is valid 
up to the second frequency moment, because in hydrody- 
namics the liquid is treated as continuum without atomic 
structure, and such an information on structure and in- 
teratomic potentials appears only within the fourth and 
higher frequency moments of S(Q,u>). By expanding the 
density autocorrelation function in a Taylor series, one 
can easily find a connection with the frequency moments 
of the dynamic structure factor as: 



d n F(Q,t) 



dt" 



t=0 



+oo 



u n S(Q,u})du 



(6) 

This relation holds for the frequency moments of both 
the collective ((uj n )s) and the self ((w")sj dynamic 
structure factor. From the previous definitions, it eas- 
ily follows that (uj°)s = S(Q) and (lj°) Ss = 1. The first 
frequency momentum are, on the other side, (lo 1 )s = 
(w 1 )5 s = ^p- and therefore are zero for any classical 
theory, characterized by symmetric spectral functions. 
The second frequency moments are (lo 2 )s — ( uj2 )s s = 
+ 0(h 2 ). Higher order spectral moments depend 
on the details of the microscopic interactions, and can be 
analyti cally derived for ad ditive pairwise interatomic po- 
tential l|De Gennesl Il959|) . Some examples will be given 
in sections III.E.5l and lTLF~2l while pratical usage of sum 
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rules for normalization purposes will be outlined in sec- 
tion nnm 



3. Quantum aspects 

The models that we will illustrate in the next sections 
have been developed for a classical system. The main 
effect of quantum-mechanical corrections stems from 
the well-known inequality of the positive and negative- 
frequency parts of the spectra, connected by the de- 
tailed balance factor e hu / k sT _ Additional sources of non- 
classical behavior, such as those associated with a finite 

1 /2 

value of the de Broglie wavelength A = (2ttH 2 /inksT) , 
are small (for lithium at melting A is only 0.11 times the 
average interparticle distance, and this ratio decreases 
for heavier metals) and can safely be neglected. Since 
the effects of the detailed balance are clearly visible in 
the experimentally measured dynamic structure factors, 
we briefly discuss a possible procedure to account for this 
quantum feature in a consistent way, while preserving the 
inherent advantages of the classical description. In doing 
this, for the sake of clarity we shall denote all the pre- 
vious classical quantities with the subscript cl, while the 
notation q will refer to the quantum case. 

The natural theoretical counterpart of the classical 
density correlation funct ion is the so called Kubo canon- 
ical relaxation function <)Kubo and Tomital . Il954() 



K g (Q,t) 



,(t) e XH 



(7) 

where /3 = l/fc^T and the angular brackets denote now 
a quantum statistical average. In the classical limit (/3 — > 
0, h — > 0) the quantum operators become classical com- 
muting dynamic al variables and K q (Q,t) — > F c i(Q,t). 
It can be shown IILovese 2 [l97l that K q (Q,t) is a real 
even function of time, so that its spectrum K q (Q, u) is an 
even function of frequency. On the other hand, the ex- 
perimental scattering cross section involves the Fourier 
transform S q (Q,u>) of the quantum density correlator 
F q (Q,t) = (l/AQEj 7 (e-'Q- ? '(°>e iQ ^ W) , The relat ion 
between S q (Q,uj) and K q (Q,oj) reads ijLovesevl Il987|) 

and, as can be easily checked, satisfies the detailed 
balance condition. Moreover, it can be seen that the 
relation 



/, ,2n\ /, ,2n— 1\ 

(U )K-p R (u> ) S 



(8) 



connects the even frequency moments of K q with the 
odd ones of S q . In addition to that, the same memory 
function framework which will be outlined in sections 



III.E.5l and lII.F.ll can be phrased for the Kubo relaxation 
function and for its Laplace transform K q (Q, s) 

By virtue of all these properties, in a situation where 
the quantum aspects not associated with detailed balance 
are marginal, it is reasonable (although not strictly rigor- 
ous) to identify the spectrum K q (Q,uj) with the classical 
quantity S c i(Q,uj) so that 



1 



S c i{Q,u) 



(9) 



Having assumed such a correspondence, from now on 
we will drop out the subscript cl and refer to the classical 
quantities as in fact done at the beginning of this section. 

The transformation © allows one to test classical 
models against experimental data. It is worth to point 
out, however, that it alters the frequency moments: as 
shown in the previous section, for instance, it introduces 
a h 2 correction to the second frequency moment, though 
this effect h as been shown to be ha rdly noticeable in liq- 
uid metals <|Scopigno et all 12000 b|) . 



B. Single particle dynamics in the hydrodynamic regime 

The time evolution of the single particle density can 
be easily obtained through the continuity equation and 
the constitutive relation (Fick's law) relating density and 
current variables: 



p s (r,t)+V- J s (r,t) = 
J,(r,t) = -DV Ps {r,t) 

It is worth to stress that while the first equation is ex- 
act, the second is a phenomenological "closure". Com- 
bining the two equations, one gets the diffusion equation 
straightforward: 

Ps(r,t) = D\7 2 Ps (r,t) (10) 
which, in the reciprocal space, has the solution 



p a (Q,t) = Ps(Q)e 



-DQ 2 t 



(11) 



The normalized autocorrelation function of the single 
particle density is now obtained as 



$ s (Q,t) = (p s (Q,t)p s (-Q)) = 



- P -DQ 2 t 



(12) 



while its fourier transform, the self dynamic structure 
factor, will read 



^'^^n^ + iDQ 2 



(13) 
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i.e. a Lorenzian function centered at lu — with 
FWHM equal to 2DQ 2 . It is worth to point out 
how, in the hydrodynamic limit, the diffusion coeffi- 
cient is related to the dynamic structure factor as D — 
limg^o fyi 7r S s (Q, u>). 

Finally, for completeness, the corresponding Van-Hove 
self correlation function is: 



G s (r,i) 



(4nDt) 3 / 2 



-r 2 /4Dt 



(14) 



p{Q,t) 
p(Q) 

S{Q,u) _ i 
S(Q) 2tt 















1 




2^ 





l' T Q 2 t _j_ }_ e -TQ 2 t 

1 

2D T Q 2 



cos c s Qt 



7 ) oj 2 + {D T Q 2 



TQ 2 



TQ 2 



(u + c s Q) 2 + (TQ 2 ) 2 (u - c s Q) 2 + (TQ 2 ) 2 

(16) 



having defined 



C. Collective dynamics in the hydrodynamic regime 

In a similar manner as in the previous section, one 
can build again a set of closed equations but, in this 
case, Fick's law does not apply, and the situation is more 
involved. The constitutive relations, indeed, couple to- 
gether the three conservation laws for the microscopic 
variables density, momentum and energy, which, in terms 
of the correspondent fluxes reads: 



p(r,t) + V- J(r,t) = 
j(r,i) + V-cr(r,f) =0 
B(r,t) + V-H(r,t) = 



(15) 



were we have defined the momentum flux o~(r, t) and the 
energy flux H(r, t) 



a a ^(r,t) = S a ^P(r, t) - T) 



du a (r,t) dup(r,t) 



dra 



dr a 



+ Mjl-OV-Kr.t) 
H(r,t) = hu(r, t) - «VT(r, t) 

here rj and £ are the shear and bulk viscosities, P and 
T are the local pressure and temperature fields, h is the 
enthalpy density and k is the thermal conductivity. 

Two additional constitutive relations (the Navier 
Stokes equation and the Fourier law) call into play the 
two additional thermodynamic variables pressure P and 
temperature T. Invoking the thermal equilibrium and 
the equations of state one obtains a closed set of equa- 
tions which can be solved to get the density density cor- 
relation function. The detailed derivation of p(Q, t) is be- 
yond the purpose of the present review, and can be easily 
retrieved in classical textbook s ijBerne and PecoraLll976t 
lHansen and McDonald Il986|) . Here we will only recall 
the final approximate results 2 which are: 



2 The exact hydrodynamics expression contains a small additional 
contribution which makes the Brillouin components asymmet - 
ric, as emphasized in iNichols and Carom el ll968tlVerkerkll200ll) . 



- 2L 
C V 

Dn 



pmCf 



2 pm 



:% + VB 



(7~1)« 



cp 



(17) 



In the above expressions cp and cy are the specific heat 
ratios at constant pressure and volume, r\ s and rjp are the 
shear and bulk viscosities, k is the thermal conductivity 

and c s — sjl i s the adiabatic sound velocity. 

The classical hydrodynamics, therefore, predicts in the 
long wavelength limit (Q — * 0) a frequency spectrum for 
the density fluctuations constituted by two main features. 
The central part of the spectrum is dominated by a quasi- 
elastic, non propagating mode related to entropy fluctu- 
ations (Rayleigh component) of linewidth T qe = 2DtQ 2 , 
which reflects the fact that thermal fluctuations decay 
over a finite lifetime r = 2/T qe . Beside, two symmet- 
rically shifted inelastic components peaked at frequency 
lo s = ±c s Q are the signature of propagating pressure 
waves (Brillouin doublet), which are damped by a combi- 
nation of viscous and thermal effects. The ratio between 
the Rayleigh and the Brillouin component is given by the 
Landau-Placzeck ratio 



Ir_ 
2I b 



= 7 - 



(18) 



Usually, the hydrodynamic regime is investigated by 
visible light scattering spectroscopy (BLS, Brillouin Light 
Scattering). In the case of liquid metals, the light scat- 
tering study of density fluctuations, is prevented by the 
non-transparent nature of these systems 3 . 



The resulting lineshape can be easily recognized as the already 
mentioned Damped Harm onic Oscillator, origina lly proposed 
within a solid-like picture and Dorneii 119971) . which can 

be actually retrieved by a liquid-like point of view within the 
memory function formalism, shown in section fll.F.2l 
In the case of copper, for instance, the distance to the Fermi 
surface is 2.3 eV. Thus, electrons are promoted by energies asso- 
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Inelastic scattering experiment can, in fact, only be 
performed by means of higher energy photons (X-rays) 
or by neutrons, but in both cases the probed wavevectors 
are fairly outside the strict hydrodynamic region. 

The way the Brillouin triplet evolves at finite Q is far 
from being fully understood, though some simplified phc- 
nomenological models h a ve be en proposed in the past 
l|McGreevv and Mitchell, 119851) . Primarily, one should 
account for the frequency dependence of the transport 
coefficients which, as we shall see, corresponds to aban- 
don the hypothesis of a Markovian dynamics. Second, 
once the wavevector approaches the inverse interparticle 
distances, structural effects are expected in the form of a 
Q dependence of all the thermodynamic quantities. Last, 
but not least, it is highly questionable whether the role 
of the thermal and viscous processes remains well sepa- 
rated at high Q. In particular, specifically in the case 
of liquid metals, due to the high thermal conductivity 
one expects DtQ 2 to become soon of the order of the 
brillouin frequency, so that entropy and density fluctua- 
tions become closely interwove . Strangely enough, this 
aspect, with a few exceptions l|Fabeit Il972j) did not re- 
ceive much attention in the past, although lately it has 
been the matter of some deb ate llScopigno and Ruoccol 
12004 ISingh and Tankeshwart 120031 I2004D . 

Before coming in the discussion of the evolution of 
S(Q,oj) at increasing Q values, above the hydrodynamic 
limit, it is worth to discuss another analytically solvable 
case: the high Q limit. 



D. The short wavelength limit 

In the previous sections we examined the self and col- 
lective motion at very small wavevectors and frequen- 
cies. In the opposite regime, i.e. at short distances and 
timescales, the particles of a fluid are expected to move 
as they were free, as it happens in an ideal gas. Since 
in this case the behaviors of different particles are un- 
correlated (Gd{r,t) = 0), the self and collective dynamic 
structure factors coincides in this limit. 

The correlation function of Eq. © can be easily cal- 
culated for a free particle, i.e. for rj(i) = rj(0) + Vjt. The 
classical free particle correlation function is: 



-iQ 



,-iQ 



(19) 



ciatcd with the blue-green end of the spectrum. As a result, red 
and orange light at the opposite end of the spectrum is reflected 
back and gives copper its characteristic color. With the alkali 
metals, the s electron is involved in promotion to the Fermi level. 
There is little overlap to the empty 3p and 3d orbitals that con- 
tribute to the conduction band. Therefore, only radiation close 
to the ultraviolet region is absorbed and visible light is reflected, 
hence the silver-like appearance of the alkali liquid-metals. 



evaluating the thermal average one easily gets 

f P 2 .„ . k B TQ 2 t 2 

F{Q,t)= j e'^^e- %Qpmt dp = e 

and the corresponding dynamic structure factor is 



(20) 



S(Q,u) 



2Trk B TQ 2 



_g 2k B l 



(21) 



By virtue of the previously mentioned considerations, 
the Van-Hove self correlation function reads: 



G s {r,t) 



2irk B Tt 2 



3/2 



(22) 



In the quantum case the correlation function (|19[) can 
be evaluated treating Vi and pi as operators, and paying 
attention to the fact that, in this case, the product of 
the exponential in Eq.© can not be reduced to a single 
exponential, as in the classical treatment. Invoking the 
identity e^e^ = e ^+^+n^\, holding when ,as in the 
present case, A, B is a number, one can write: 



e -iQT 4 (0) e iQT 4 (t) 



«»•?{■*-¥[?!■*.'«(*)] 



ink 



and, using the result of Ea. (|19|l . one obtains: 

F(Q,t) = e ~^ kBTt2 - lht ) (23) 
the correspondent dynamic structure factor: 



S(Q,co) 



2irk B TQ 2 



2k B TQ-* 



- 



(24) 



Summing up, the quantum dynamic structure factor 
for a free moving particle is a gaussian with recoil energy 

wr(Q) = ^ and linewidth a = \p^f-Q- It is worth to 
point out that Eo.l)24f satisfies the detailed balance con- 

dition S(Q,uj) — e k B T S(Q, — uS) and coincides with the 
classical case for h — > or, equivalently, for T — * oo. As 
far as the sound velocity is concerned, one can still define 
the apparent frequencies u>i(Q) = ^(ujr± a/8<t 2 + uij^) as 
the positive and negative maxima of the longitudinal cur- 
rent Cl(Q,ui) = u S gi' u ' . In the classical case, the two 
values coincides and are <jJi{Q) = y/2a: 
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E. The non-hydrodynamic region: single particle 



1. The gaussian approximation 

As can be easily noticed looking at the expressions l|14|) 
and 11' 11 , both the hydrodynamic and the ideal gas limit 
end up with the Van Hove correlation functions that are 
gaussian in r. On the basis of this observation seems 
natural to assume the gaussian dependence as valid in the 
whole dynamical range. In terms of the second moment 
of G s (r,t) one can write the following expression: 



G s (r, t) 



2<r»(t)) 



2(r^(t)> 



(25) 



where (r 2 (t)) is the mean square displacement which, 
in the hydrodynamic and single particle approximations 
reads (r 2 (t)) = 6Dt and (r 2 (t)) = ^ft 2 , respectively. 

In the gaussian approximation, therefore, the self scat- 
tering function is related to the mean square displace- 
ment which can be, for instance, inferred by molecular 
dynamics simulations. 



2. The jump diffusion model 

The jump diffusi on model was firstly intr oduced by 
Chudley and Elliot l|Chuldev and ElliotlTl96l|) . The par- 
ticle is thought to live for a residence time To in the cage of 
its neighbors, and at some point to change cage. In some 
sense, therefore, it is the opposite of collisional models, 
where the free diffusion of a particle is sometimes inter- 
rupted by collisional events. The jump diffusion model 
sets a rate equation for the Van-Hove self scattering func- 
tion of the kind: 



5G s (r,<) 1 



dt 



Afro 



with Af the number of available residence sites. 

By fourier transform in space and time one immedi- 
ately gets: 



Q(n , i HQ) 



nu 2 + f{Qf 



(26) 



which is a lorenzian function with Q— dependent damp- 
ing 



f(Q) = -— E^' 1 - 1 ) 



TLTQ 



(27) 



This latter can be conveniently estimated supposing that 
the vectors 1 have random and continuous orientations 
and distributions. In this case one can average Eq. l|27|l : 



f(Q) 



i 

TO 



1 - 



1 + Q 2 l 2 



It can be easily shown that Eq.l| 2 ^[l tends to the Fick's 
free diffusion expression JT3J (Egclstafi[llH3) 



3. The mode coupling theory 

The Laplace transform of the intermediate scattering 
function can be generally written as: 



Ss(Q,s) 



} 2 U{Q,s) 



being U{Q,s) a generalized frequency and wavevector de- 
pendent diffusion coefficient. The mode coupling the- 
ory provides a self consisten t expression for U s (Q,s) 
l(de Schenoer and Ernsti li"979T) . and the resulting self dy- 
namic structure factor reads: 



S s (Q,u) w - 



DO 2 



1 



ttuj 2 + (DQ 2 ) 2 ttDQQ 



-ReG 



ilo + DQ 1 



SDQ 2 



with 



G(s) = arctan 



(s - 2)V~s 



(28) 



(29) 



where Q* = WnmpD 2 /k B T and 6 = D/(D + u), being 
D the diffusion coefficient and v the kinematic viscosity. 

An estimate of the FWHM can be numerically evalu- 
ated l)De JoneLll99^1 . yielding: 



Q 



H{5) + Q(Q 3/2 ) 



DQ 2 



(30) 



with 



H(5) ra 1.45<5 3/2 [1 - 0.735 - 0.15<5 2 - Q{5 3 )] (31) 



4. The Nelkin-Ghatak model 

Nelkin and Ghatak have considered a dilute gas in 
which the atomic motion is dominated by binary colli- 
sions, with a distribution function obeying a linearized 
Boltzman's equation, valid in a small disturbance limit, 
i.e. for arbitrary large fluctuations compared to the 
mean collision time. In terms of the reduced variables 
x = —uj/Qvo and y — a/Qvo, with u — y/2KT/m and 
a an adjustable parameter, and by introducing the real 
(u(x,y)) and imaginary (v(x,y)) parts of the probability 



12 



integral for complex argument z(x + iy) 
t^dt 



e (z ■ 



U(x,y) 
V(x,y) 



Vnyu(x, y) 
y/%yv(x,y) 



(32) 
(33) 



one gets the foll owing expression for the incoherent 
scattering function ijNelkin and Ghataktll964|) : 



1 U{l-U)-V 2 
na (1 - U 2 ) + V 2 



(34) 



it can be easily noticed that Eq. l|34|l has the correct low 
Q (lorenzian) and high Q (gaussian) limits: in the first 
case it is sufficient to pose a = v 2 /2D, while at high Q's 



one has that y 
recovered. 



and the familiar gaussian shape is 



being — due to the orthogonality of p s and p s . 

A good description of the evolution of cj) s {Q 1 1) is nor- 
mally gained utilizing the first two equations of the chain 
QH5JI. In terms of Laplace trasform it holds: 



4> s (Q,t) = f« + P(Q,«) 



s + 



Af«(Q,t = 0) 



s + M(2)(Q,s) 



(38) 

The initial values M^(t = 0), can be easily obtained 
from the general relation: 



M w (0) = -M^-^O)- [M (i_1) (0) 



(39) 



which, in turns, is obtained deriving Eq.||35j) and exploit- 
ing Ea. l|36|) . For the first two memory function it holds 
therefore: 



5. The memory function formalism 

The easiest way to abandon the hydrodynamic region 
is to assume the frequency dependence of the transport 
coefficients, entering the so called generalized hydrody- 
namics. The natural playground for performing such step 
is the memory function framework IjMorl 119651) : we will 
recall here the basic formalism while a detailed treatment 
can be found in specialized b ooks (Balu cani and ZoppI 
Il983t lHansen and McDonald I1986J) - Let be the 

correlation matrix of a given set of dynamical variables A 
(M„°J = (A* v A a (t))). The equation of motion of M(°>(f) 
can be conveniently expressed in terms of a chain of ar- 
bitrary order n of integro-differential equation involving 
appropriate memory functions M'''(f) for i — l...n: 



M W(Q,o)^ = (- 2 W)> ; 



(Q)h 



(40) 



Here (tu n (Q))s s are the frequency moments of S s (Q,uj) 
and the quantity fig is related to the mean squared force 
(|F| 2 ) acting on the diffusing particle and, for a system 
of identical particle interacting via pairwise interactions 
potential V(r), it holds: 



P 

3m 



V 2 V{r)g{r)dr (41) 



dt 
with 



Jo 



(35) 



(36) 



here fi^ -1 ) is a set of generalized frequency matrixes, 
while the memory kernels M l (i) rule the dynamical evo- 
lution of the observables correlation matrix M^ ^. 

In the specific case of self dynamics, as we have seen 
in section (|II.BI| the relevant set of variable is given by 
the self density only therefore: M^(t) = <j) s (t). The 
equation of motion for the density correlation function 
will be: 



d^ s {Q,t) 
dt 



+ M^(Q,t-t')(f> s (Q,t')dt' = (37) 



where g(r) is the pair distribution function. 

It is worth to stress that Eq. (|38|l is the exact solution 
of motion in which all the dynamics is detailed by the 
shape of the second order memory function M^(Q,t). 
The most common way of solving the equation is making 
a guess on the shape of the memory function. A use- 
ful approximation is provided by the exponential shape, 
which has the advantage of being easily Laplace trans- 
formed: 



M^(Q,t) = A 2 (Q) e - t ^W) = [2(u 2 (Q)) Sa + O 2 ] e"*/^« 

(42) 

With such a choice it follows straightforward: 



S s (Q,u) = - 



( W 2 ) Ss (2<c 2 ) Ss +^)r s 



7TLJ 2 t s (lj 2 - 3(io 2 ) Ss - n 2 ) 2 + (uj 2 - (cu 2 ) Sa ) 2 

(43) 

It is interesting now to look at the FWHM T S (Q) of 
Ea. (|43|) . In the small Q limit it is easy to show that 
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r,(Q) 



1 



1 



2knT Q 2 



lil'r, 



(44) 



The memory function approach with exponential ker- 
nel, therefore, predicts a quasielastic lineshape narrower 
then the hydrodynamic one, a result which is in agree- 
ment with several experimental data. Contrarily, in the 
gaussian approximation the linewidth is always larger 
than the hydrodynamic value. 

Theoretical expressions have been proposed in the 
past for the memory function in terms of kinetic the- 
ory, splitting the memory function in a contribution due 
to uncorrelated binary collision, obtained by a Fokker- 
Plank equation, and a long time contribution repre- 
senting the coupling of a tagged particle to the col- 
lective motion of the surrounding particles. This ap- 
proach have been tested in hard spheres, Lennard-Jones 
and alkali metals llBengtzelius et a/lll984l:ISiogrenlll979t 
Sjo gren and Siolandeii Il979|) 



F. The non-hydrodynamic region: collective motion 

1. The Langevin equation 

The most reliable approach to the study of collective 
dynamics at finite wavevectors parallels the one adopted 
for the single particle in section III. El i.e. an extension 
of the classical hydrodynamics assisted by the formalism 
of the memory function ruling the Langevin equation of 
motion of the density fluctuations. 

We will deal, therefore, with a 3 x 3 correlation matrix 
M'W, and a set of n (with n arbitrary large) memory 
matrixes with the same dimensionality and cou- 

pled by the chain of n equations (|35|l . Actually, in or- 
der to work with a orthogonal set of variables (and en- 
ergy and density are not), one normally prefers to replace 
the microscopic energy with the microscopic temperature 



T{Q) 



E(Q) 



(E'(Q)p(Q)) (n] 

( P '{Q)p{Q)) p[ - Li ' 



With this 



mpc v (Q) 

choice, solving for cf>(Q,t) 
F(Q,t)/S(Q) in terms of Laplace transform one has: 



Mj$(Q,t)/M$(Q,Q) = 



s) = 



KQ) 



s+M { j)(Q,s)- 



(MY'(Q.s)- 



where 



= k B TQ 2 /mS{Q) 



(45) 



(46) 



seems natural to introduce a finite Q generalization of 
the isothermal sound velocity as c t (Q) = uiq(Q)/Q. 

It is easy to recognize in Eq. (|45|l the same structure 
of Ea. l|38|l . once the following identification is made: 

_ (M^(Q, a) - iSl%(Q))(M$ (Q, s) - iSl ( °l(Q)) 

The effective memory function Mpp (Q, t), therefore, 
can be formally considered as the second order memory 
function of a chain of two equations for the single variable 

P (Q,ty. 



+ / M^(Q,t-t')F(Q,t')dt' = (47) 



dF(Q,t) 
dt 

dM(jQ,t) + f' M (2) (Qjt _ /)M ( 1)(Q , t ')dt' = 
dt Jo 

which, as can be easily verified, corresponds to the single 
second order integro-differential equation: 



</>(Q,t) +w o 2 (O)0(Q,t) 



+ / M(Q,t-t') 4> {Q,t')dt' = 



(48) 



From here on, to save writing, we define M(Q,t) = 
M^(Q,t) = M^ ff) (Q,t).^ 

From the knowledge of <fi(Q, s) one straightforwardly 
obtains S{Q,uj) = {S (Q) / ir}dt{(f)(Q , s — iuj)} in terms of 
the real (M 1 ) and imaginary (— M") parts of the Fourier- 
Laplace transform of the memory function: 



S(Q) 



w§(Q)M'(Q,w) 



S(Q,u) 

* [lo 2 - w§ - uM"(Q, uj)Y + [wAP (Q, uj)Y 

The spectral features of the dynamic structure fac- 
tor can be characterized by its frequency moments 
(uj n (Q))s = f uJ n S(Q,uj)did, where, for a classical 
system, only the even frequency moments (such as 
(u°(Q))s = S(Q) and (uj 2 (Q)) s = Q 2 are different 
from zero. 

It can be easily proven that the dynamic structure 
factor is related to the longitudinal current spectrum 
through the relation C L (Q,u) = (uj 2 /Q 2 ) S{Q,u). The 
presence of the factor uj 2 wipes out the low frequency por- 
tion of the dynamic structure factor, and consequently 
emphasizes the genuine inelastic features of S(Q,u>). Af- 
ter its definition and Ea. (|38|l . it is readily seen that the 
Laplace transform Cl(Q, s) satisfies 



Recalli ng the defi nition of the isothermal sound veloc- 
ity c t = y/l/pmXfr, and the expression for the low Q — > 
limit of the static structure factor S(0) — pxt^-bT it 



C L (Q,s) 



-s[sF(Q,s 
k B T 



S(Q)} 



(50) 



{ S + [lo 2 (Q)/ s } 



M(Q,s)} 1 
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Again, the spectrum Cl(Q,oj) can be expressed as 
(1/7t)3?{Cl((5, s — iw)}. Then the position and the 
width of the inelastic peaks in Cl(Q,uj) are determined 
by the poles of Cl{Q, s). 

2. Collective memory function and hydrodynamics 

The effective, second order memory function M(Q,t) 
accounts for all the relaxation mechanisms affecting col- 
lective dynamics and, consequently, is the central quan- 
tity in most theoretical approaches. In analogy with the 
single particle case (Eq. (00))), the initial value of M(Q, t) 
is related to the spectral moments of S(Q,ui) by: 



M(Q,0) 



( W 2 (Q)) S = A 2 (Q) (51) 



Along the same line, relations similar to l|40ll and 141|) 

holds, and an explicit expression for (lu 4 (Q))s can be 
given, involving both the derivatives of the interparticle 
potential and the pair distribution function: 

(^ 4 (Q)) S 3K B TQ 2 , p f d 2 V(r) 



2 (Q)K< 



dz 2 



-(l-e- lQz )g(r)d 3 r 

The second and fourth frequency moments are par- 
ticularly significant, as they rule the sound velocity in 
the whole Q — lj domain. For sufficiently large s, in- 
deed, M(Q,s) ~ M(Q,t = 0)/s and Eq. JSDJ} is seen 
to have poles at s = ±i^uj 2 (Q) + A 2 (Q) = ±iuj OQ (Q). 
This latter relation defines the frequency which 
characterizes the instantaneous collective response of the 
liquid at the wavevector Q and, in turns, defines the un- 
relaxed sound velocity as c<jo(Q) = lu qo (Q)/Q. In the 
opposite limit, i.e. for small s, one easily verifies that 
M(Q,s) ~ J M(Q,t)dt. In this (relaxed) regime, there- 
fore, the poles of the longitudinal current are located at 
s = ±iwg(Q), i.e. the longitudinal modes propagate with 
the isothermal sound velocity ct(Q) — u)q(Q)/Q. Sum- 
ming up, whatever the details of the memory function 
are, in presence of a relaxation process one observes a 
transition of the longitudinal sound velocity between two 
different regime, associated to the evolution of the lon- 
gitudinal current correlation maxima. It is interesting 
to study in a parallel way the evolution of S(Q,u>). In 
the low frequency limit, Eq. 03 reduces to the spectrum 
of a damped harmonic oscillator (DHO) of characteristic 
frequency ujq (which, in general, does not coincides with 
the position of the inelastic maximum), with damping 
r = / M(Q,t)dt. 



with an higher characteristic frequency ±o; 00 (Q). Addi- 
tionally, it appears an elastic peak of area 



V(Q) 



wg(Q)+A=(Q)- 



The same results can be directly retrieved in the time 
domain from Eq. I|49(l . substituting the high frequency 
limit (constant) and low frequency limit (delta shaped) 
of the memory function, in turn, and taking the Fourier 
transform. 

The memory function formalism leads, in the small 
wavevectors limit, to the hydrodynamic prediction Of 
Ea. ((TT)|) . It can be easily proven, indeed, that the diag- 
onal terms of the memory matrix appearing in Eq. i|45|) 
have a Q 2 dependence, while the cross terms follow a 
Q 4 dependence. For Q — > 0, therefore, these latter can 
be neglected. Moreover, from the continuity equations 
(|15|l it follows that in the same limit the time depen- 
dence of the conserved quantities and their associated 
currents becomes increasingly slow. Consequently, one 
can model the decay of the terms Mjj(Q — > 0,i) and 
M^t(Q ~ * Oj t) as instantaneous, or equivalently, con- 
stant in the Laplace domain. By making the identifi- 
cations M}j{Q -> 0, s) = D V Q 2 and M\ T (Q -> 0, s) = 
jD T Q 2 7 and computing n^(Q)fi^(Q) = ( 7 -l)wg(Q), 
from Eq. (|45|l one gets for the memory function: 



M(Q^0,s) = D V Q 
M{Q -> 0,i) 



(7-iH 2 (Q) 



s + 1 D T Q 2 
2D v Q 2 5{t)+Lol{Q){ 1 - 



(54) 



The dynamic structure factor can be obtained by sub- 
stituting Eq. (|54|l in the general expression Ij49|l . In 
the Q — ► limit, moreover, one has loq — ctQ >> 
-yD T Q 2 . In this limit, S(Q,w) is: i) a DHO func- 
tion around the Brillouin peaks, which are located at 
lu s = \p^c t Q (adiabatic sound propagation) and have a 
linewidth DyQ 2 + (7 — 1)DtQ 2 ii) a lorenzian function 
around to = 0, whose linewidth is 2DtQ 2 - In the small 
damping limit, the DHO is well approximated by two 
symmetrically shifted lorenzian, and the hydrodynamic 
limit of Ea. 1)16(1 is finally recovered. 

The advantage of the memory function approach is, 
however, in providing a way to generalize the hydro- 
dynamic result for wavevector and frequency dependent 
transport coefficients. To this purpose, from the very 
start it is convenient to separate in M(Q,t) the decay 
channels which explicitly involve couplings to thermal 
fluctuations ( M t h(Q,t) ) from those directly associated 
with longitudinal density modes ( Mi,(Q,t) ). 



S(Q,u) _ 1 co 2 (Q)T(Q) 
S(Q) it {uj 2 - wg) + uj 2 T 2 



(53) 



In the opposite, high frequency, limit, again Eq. I|49|) 
tends to an harmonic oscillator with no damping and 



3. Finite wavelengths generalization: beyond hydrodynamics 

A straightforward generalization of ordinary hydrody- 
namics at finite wavevectors suggests for the thermal con- 
tribution the following form 
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M th (Q, t) » ( 7 (Q) - l)wg(Q) exp[- 7 Z? T (Q)Q 2 t] (55) 

where Dt(Q) and 7(Q) can be regarded as a finite Q 
generalization of the quantities Dt = n/nCp, being k 
the thermal conductivity, and 7 = cp/cy. 

It must be stressed, however, how the extension to fi- 
nite wavevectors, in the special case of highly conduc- 
tive systems, requires attention to the physics b ehind this 
model. As pointed out by Faber ijFabert Il972jh in a liq- 
uid metal, owing to the high thermal conductivity, the 
quantity jD T (Q)Q 2 may easily become larger than the 
Brillouin frequency as soon as Q » Dt^Q) 4 ' Actually, 
both Ct(Q) and Dt(Q) are expected to decrease on ap- 
proaching values comparable to the inverse mean inter- 
particle distance, i.e. in coincidence with the first sharp 
rising edge of S(Q). But assuming that the transition 
occurs below this edge, one finds a crossover condition 
of Q — -jf- which, considering typical values of sound 
speed and thermal diffusivity of metals (a few thousands 
meters/second and ~ 50 nm _2 /ps, respectively), lies at 
wavevectors around 0.1 4- 0.5 nm _1 , which is indeed con- 
sistent with the initial assumption. In other words, on 
increasing Q's, the thermal peak broadens ultimately 
overlapping with the Brillouin lines, the sound propaga- 
tion turns from adiabatic to isothermal, and independent 
thermal fluctuations become impossible. By the mem- 
ory function point of view, at wavevectors Q » 
(i.e. when the condition LOT t h « 1 holds), M t h de- 
cays instantaneously in a similar fashion to Ml- In this 
condition, the dynamic structure factor is a DHO with 
isothermal characteristic frequency, while the brillouin 
damping is given by the area under the memory function 

which is DyQ 2 + 7 rj^ ■ The existence of this adiabatic 
to isothermal crossover, is expected just below the lower 
accessible Q's of an inelastic scattering experiment (« 1 
nm _1 at present) and therefore has never been observed 
directly. Moreover, the magnitude of this effect is ruled 
by 7 — 1, which is normally very small. A remarkable ex- 
ception to this latter condition is constituted by Nickel 
(Bcrmeio et alll200Cft (7 = 1.88), which has bee n studied 
by INS, as it will be discussed in the next section. 

Beside the thermal contribution, also the viscous term 
Dy is expected to exhibit a Q dependence. For simpli- 
fied Lennard-Jones pairwise interaction it has been shown 
a decay of rji{Q) of a factor ten up to the main peak 
of S(Q), with 77b going negative over the same region 
(|TankeshwarMl994j) . 



4. Finite frequencies generalization: viscoelasticity 

A second important generalization of transport coeffi- 
cients concerns their possible frequency dependence. 

This latter case stems when the frequencies of the ob- 
served density fluctuations are high enough that their 
timescale competes with the one ruling the decay of Ml,. 
In this condition, one has to drop the hypothesis on the 
instantaneous (Markovian) nature of the viscous term, 
and has to introduce a finite timescale (r) for the decay 
of Ml . As discussed at the beginning of this section, the 
timescale of Ml sets a new crossover between two differ- 
ent regimes, characterized by different sound velocities 
and attenuations. 

The simplest practical way to go beyond the hydrody- 
namic result l|54l) is to allow for an exponential decay of 
M L (Q,t): 

M L (Q,t) = A|(Q)e-*/^) (56) 

with A|(Q) = lol{Q) 2 — J^oiQ) 2 , in order to have the 
correct normalization of the whole M(Q,t = 0). 

Although this has the advantage of analytical simplic- 
ity when dealing with Fourier transform, a drawback of 
this ansatz lies in the violation of some basic short-time 
features of the memory function (such as a zero derivative 
at t = 0), causing the divergency of (ui n )s for n > 6. 

Neglecting thermal effects, Eq. l|5T)|) yields the so - 
called viscoelastic model 5 for S(Q.lj) l|Lovesevl 1 1 9 7 lh . 
Since as Q — > Ml(Q,s = 0)/Q 2 can written as 
[c 2 ^ — Cq]t(Q — > 0), the requirement that this coincides 
with r]L/nm shows that the time t(Q) must be finite as 
Q — ► 0. Such a connection with viscous effects justifies 
the physical interpretation of the rate l/r(Q) as a param- 
eter giving an overall account of all relaxation processes 
by which the longitudinal response of the liquid is af- 
fected by time-dependent disturbances. In particular, for 
slow perturbations developing over a timescale t 3> t(Q) 
the system can adjust itself to attain local equilibrium 
and the usual viscous behavior holds. In contrast, for 
density fluctuations fast enough that t -C t(Q) the liq- 
uid responds instantaneously, with a solid-like (elastic) 
behavior. The crossover between these limiting situa- 
tions marked by frequencies to such that lut(Q) sa 1) 
is ultimately responsible for the gradual changes often 
detected in the sound dispersion of several liquids at in- 
creasing wavevectors. 

Although appealing, the simplicity of the viscoelastic 
model can be deceptive. First of all, the model itself pro- 
vides no clue for the physical origin of the decay mecha- 
nisms leading to the rate 1/t(Q). 



4 More specifically, the discrepancy observed in liquid lead be- 
tween the sound velocity measured with ultrasound (adiabatic) 
and with inelastic neutron scattering has been tentatively as- 
signed to the isothermal nature of the sound propagation at the 
wavevectors probed with neutrons. 



Actually, the viscoelastic model stems from the approximation 
. with the additional condition 7 = 1 (Afj^(Q,f) = 0). 
Within the viscoelastic framework, indeed thermal effects are ne- 
glected, in the sense that the hydrodynamic limit is isothermal 
(i.e. Al(Q) = co L (Q) 2 -w (Qf) 
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Actually, the situation is even more involved. In ear- 
lier MD studies of Lennard Jones fluids, it was soon real- 
ized that the viscous dynamics in the microscopic regime 
(i.e. at wavelength comparable with the inverse mean 
inter-particle separation) proceeds through two distinct 
processes, characterize d by two well separate time scales 
l|Levesaue et al\Tl97?§ . More recently, the advent of IXS 
provided the experimental evidence substantiating these 
speculations jMonaco et a71 . l2004l:IScopigno et ad 12000c. 
200ll l2002albft . 

In view of these result, the obvious remedy is to modify 
the simple ansatz (|56|) by allowing a two step decay of 
M L (Q,t): 



M L (Q, t) = A|(Q) (1 - a(Q))e^ (Q)t +a(Q)e-^ Q)t 



where the rate 71 (Q) is chosen to be larger than 72 (Q), 
so that the dimensionless factor a(Q) measures the rel- 
ative weight of the "slow" decay channel. Besides being 
more flexible than the viscoclastic model, we shall see 
that the ansatz l|57|) has the much more important merit 
that the presence of two different timescale does in fact 
have a definite physical interpretation. 

A simplified version of the previous ansatz has been im- 
plicitly introduced in the viscoelastic analysis of Brillouin 
Light Scattering (BLS) spectra o f glass forming materials 
(see for example l)Li et all \1992). In fact, in these works, 
the general expression of Mij(Q,t) for a two times decay 
is always expressed as 

M L (Q,t) = 2Al(Q)r^Q)5(t)+Al(Q)e- t ^^ 

with explicit reference to the so-called a- (structural) 
relaxation process as responsible for the long lasting tail, 
and to the fx— microscopic process as additional, faster, 
relaxation dominant over a very short timescale (in the 
BLS window the condition lut^ << 1 holds). 

On the contrary, the above mentioned IXS works it 
have shown how this approximation is no longer tenable 
in the case of liquid metals at the IXS frequencies. In 
particular, it has been pointed out that the slower (-a) 
relaxation time satisfies the condition ujb{Q)t~u{Q) » 1, 
i.e. some part of the viscous flow is frozen. As a conse- 
quence, at the wavevectors typical of the IXS experiments 
(Q = 1 -7- 20 nm -1 ) the quasielastic spectrum acquires a 
component arising from this frozen structural relaxation. 

The origin, at the atomic level, of this fast decay chan- 
nel is still an open issue: the rapidly decaying portion 
of Ml(Q, t) is customarily attributed to largely uncorre- 
lated collisional events, similar to those occurring in a di- 
lute fluid. In addition, at the high densities typical of the 
liquid state, non-negligible correlations among the colli- 
sions can be expected, making no longer valid an inter- 
pretation only in terms of "binary" collisions. Although 
the magnitude of the correlation effects is relatively small 
and their buildup slow, once established their decay is 
even slower, and for t > l/7i(Q) = t (Q) this relax- 
ation channel dominates the decay of Ml(Q,£), which 



consequently may exhibit a small but long-lasting " tail" . 
The ansatz (|57f) can incorporate most of this physics: 
on the basis of the latter, one may reasonably antici- 
pate that a(Q) -Cl, and that the time l/j2(Q) = T a {Q) 
is distinctly longer than l/7i(Q) = T fl (Q). In this pic- 
ture, the best fitted values of the viscoelastic rate l/r(Q) 
clearly represent some sort of "weighted average" be- 
tween 71 (Q) and 72 (Q). Finally, we may argue that at 
increasing Q (namely, over a shrinking length scale) the 
magnitude a(Q) of correlation effects should decrease, 
and that at higher temperatures the value of a(Q) at a 
given wavevector should equally decrease. On a general 
basis, the requirement that limg^o Ml(Q, s — Q)/Q 2 — > 
rjL/nm now takes the form 



(4 - eg) 



(l-a(Q-fO)) , a(Q^O) 



7i(Q^0) 



72(Q^0) 



r} L /nm 
(58) 

The refined model l|57|l yields a dynamic structure fac- 
tor given by 



S(Q,u) = 



xRe 



S(Q) 



iw+7i(Q) 



A 2 JQ) , 
iw+72(Q) iui+ 1 (Q)D T (Q)Q' A 



59) 



One of the major drawbacks of expression (|59|l . is that 
one normally overestimates the relaxation strength of the 
faster viscous process. Such problem is somehow ex- 
pected as one is trying to force an exponential depen- 
dence to reproduce the memory function at short times, 
which, instead, has a zero derivative in the t — > limit. 
As a consequence, adjusting the characteristic time of 
the exponential memory on the experimental data one 
can reproduce the decay of the true memory function 
but will overestimate the short time limit, due to the 
cusp behavior of the exponential at t = 0. 

An obvious remedy it is a better choice of the 
memory function model. An alternative possibility 
could be a gaussian shape. Also the sech function 
has been proposed as solution of the Mori Equation 
ijTankeshwar and PathakL Il994(l . but this latter case is 
of limited practical interest, as the fourier transform is 
related to the digamma function and, therefore, the ex- 
pression H59|) must be numerically evaluated. 



G. Kinetic theories: the hard sphere approximation 

Special attention has been devoted in the past 
to the theoretical and numerical study of the hard 
sphere model, as it conveniently mimics the behavior 
of more realistic sim ple liquids ijFurtado et all Il975l 
iLebowitz et al\.\ l969\\. In the eighthes, transport coeffi- 
cients (lAllev and Aldeilll983ft and neutron scattering rc- 
sponse l|Allev*^r^/TTT98^ ) have been evaluated by means 
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of molecular dynamics. On the theoretical side, the dy- 
namical properties have been investigated by a revisiting 
the so called "Enskog fluid", i.e. by mean s of the spec- 



1984 




1984 





dc Schepper and Cohen, 1980: dc Schepper et al. 
Kamgar-Parsi et all Il987t iMrvglod et all Il995h 



This latter approach turned out to b e particularly usefu l 
to describe INS experimental data ijCohen et al\ . \l98l\i . 
and we will briefly recall the basics in this section. 

The main idea beyond the Enskog's theory is to eval- 
uate the correlation functions J2J replacing the Liouville 
operator C and the dynamical variables a a (Q), defined at 
the N— particle ensemble level, with the one particle En- 
skog's operator L, and appropriate one particle variables 
4>a(Q)- It can be easily recognized, indeed, that 



where C(Q) = [1 - ^y] and C (Q) = lim p ^ C(Q). 

An explicit expression for the dynamic structure factor 
can be easily retrieved thought the spectral decomposi- 
tion of L: 



L(Q,vi) = -5^|* J -(Q J vi)z i (QK$ i (Q J vi) 



in which Zj, and $j are eigenvalues, left and right 
eigenfunctions of — L, respectively. S(Q,lu) then reads 



ttS(Q, lj) 
S{Q) 



Re 



= Re J2 



iuj — L(Q, vi) / 1 
Bj(Q) 



iu + z 3 (Q) 



(63) 



with 



where, for the first three dynamical variables it holds 



S j (Q) = (^(Q,v 1 )) 1 ($*(Q,v 1 )) 1 



(64) 



(j>2 (v) 



S{Q) 



m Q • v 

k B T Q 
3 — mv 2 /ksT 
76 



For an hard sphere s ystem, a possible asymme tric rep- 
resentation of L reads i|de Schepper et a/1 11984^1 : 



HQ, vi) = «Q • vx + pg(*)A(Q) + P A(Q) (60) 

where g{cr) is the pair distribution function evaluated at 
the contact point between two spheres. The first term 
appearing in (|60|) is a free streaming contribution. The 
second term accounts for binary collisions through the 
operator A(Q), defined through its action over a generic 
function of the velocity /(vi) 



A(Q)/(vi) = -a 2 J da J dw^{v 2 )56{8) 

{/ (vi) - /(vi) + e-^[f(v 2 ) - /K)]} (61) 

with £(v) the normalized Maxwell distribution function, 
9{x) the Heaviside step function, a the unit vector, S = 
(vi— V2)-(T and v[ 2 = vi^T^^ the post-collision velocity. 
The third term, finally, is a mean field operator defined 
through 

M(Q)/(vi) = [C(Q) - g(a)C (Q)} 

dv 2 !i(v 2 )iQ ■ v 2 /(v 2 ) (62) 



The subscript explicitly indicates that we are now 
dealing with single particle averages over the Maxwell 
Bolzmann velocity distribution function. 

There are several approaches to determine the spec- 
trum of L, but the main point is that different approxi- 
mations can be performed according to the density and 
kinematic region of interest. These regions are marked 
by the values of the reduced density Vq/V = pa 3 /^/2, Vo 
being the closed packed volume for spheres or radius er, 
and the Enskog mean free path l& = lo/x with Iq the 
Bolzman mean free path Iq = \j \f2p-Ka 1 and \ ~ 9( a ) 
the pair distribution function evaluated at the contact 
point between two spheres. 

The lower three eigenvalues of L always go to zero with 
Q — > 0. In the same limit it can be shown that these 
latter eigenvalues are: 



zi(Q) = zh(Q) — DteQ 2 

Z2,a(Q) = z±(Q)=±ic Q + T E Q 2 (65) 

and only the first three coefficients (|64|l are relevant, so 
that the eigenfunctions are linear combinations of the 
density, current and energy variables. In other words the 
hydrodynamic result of Eq. I|16l) is recovered, with the 
dynamic structure factor composed of three Lorenzian 
functions: a diffusive heat mode and two propagating 
modes with the adiabatic sound velocity, with dampings 
Dte and Te as given within the Enskog's transport the- 
ory. The relative intensities of the thermal and acoustic 
contributions are ruled by the Landau Placzek ratio. 

This limit is practically attained at low densities 
(Vq/V < 0.1 and therefore Ie ~ ^o) an d sufficiently 
small Q's (Qa << 1), when the contribution (|62|l can be 
safely neglected and the term (|61|l can be replaced with 
A(Qa = 0). This normally occurs up to the case of light 
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scattering of dilute gases (QIq « 1), and one speaks in 
terms of tree extended hydrodynamic modes. Still in the 
same low density-momentum range, but at wavevectors 
1 < QIq < 3 a description in terms of a few hydrody- 
namic modes is no longer allowed, while for QIq > 3 one 
can evaluate the first order corrections to the free particle 
result of Ea. (|21[l which reads: 



S(Q,w) = 



'^^+ y Z +O(l/Ql 
Qh 

(66) 

with S B (tu/Q) the leading correction to the free stream- 
ing term — iQ ■ vi of Eq. I|6U|) due to a single bi- 
nary collision event, which can be numerically estimated 
l|Kamgar-Parsi et adll987ft . 

Conversely, for dense fluids (Vo/V > 0.35), the hydro- 
dynamic scheme breaks down at as soon as QIe > 0.05. 
Above this value, only the extended heat mode is well 
separated from all the other modes. Still a description 
in terms of three effective hydrodynamic modes apply, 
and the low Q limit of these modes is again coincident 
with the hydrodynamic result. For < QIe < 1 the 
free streaming and the mean field contributions of the 
Enskog operator can be treated as perturbation to the 
binary collision term (|61|l . and the following approximate 
expression for the extended heat mode can be given: 



z h (Q) 



DeQ' 
S(Q) 



d(Q) 



(67) 



in which 
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As in the low density case, for QIe > 3 the the 
free streaming limit is recovered along with the binary 
collisions corrections. In this case, however, a differ- 
ent approximation holds, since the binary collision term 
,2pn} is now replaced by the Lorentz-Boltzmann operator 
h(Q — > oo). One ultimately gets an expression similar 
to in which S b {uj/Q) is replaced by the different 
Lorentz-Boltzmann expression. 

At intermediate densities, finally, again the hydrody- 
namic result does not hold for QIe > 0.05. Moreover, in 
this regime three effective modes do not suffice and one 
has to extend the description with two additional kinetic 
modes, i.e. including the kinetic part of the z-z com- 
ponent of the stress tensor and the z -component of the 
heat flux ijKamgar-Parsi et a/.L HQS? 1 ). Like in the high 
density case, for 1 < QIe < 3 all the modes are closely 
interwoven, while at larger Q's the free streaming limit 
with the Lorentz-Boltzmann correction is retrieved. 

According to kinetic theory, therefore, S(Q,oj) in sim- 
ple liquids not too far from the melting temperature (such 
as Argon), should be described in terms of three loren- 
zians up to relatively large wavevectors (Q 30nm _1 ). 
Consequently, sound modes should exist even in a Q re- 
gion where side peaks are not distinctly observable. In 
this respects, Lovesey argued that the extended hydrody- 
namic picture should break up above Q « 3 nm " 1 , while 
abov e a viscoelastic theory should be utilized ijLovesevt 
I1984D . 

Exploiting the previous results obtained for the coher- 
ent case, one can describe the incoherent dynamics via 
the Lorenz-Enskog operator 



is the Enskog diffusion coefficient and 



L 8 (Q, vi) = iQ ■ vi + pg(a)A(Q -> oo) 



(69) 
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(68) 



can be approximated in terms of the first two even 
spherical Bessel functions. Enskog's diffusion coefficient 
is related to the Bolzman diffusion coefficient 



Da 



8pa< 



k B T 0.216 k B T 



pa* 



by the collision enhancing term g(u) as D B — Do/g(a). 
Using the analytic expression of g{cr) for non attractive 
hard spheres one finally gets an expression for De in 
terms of the packing fraction ip — irpa 3 /6: 



As already mentioned, at large Q's the self and the 
collective dynamics coincide, and therefore the spectrum 
of L will tend to the one of L s (A(Q — > oo) = 0). In this 
limit the extended heat mode Zh will tend to the self- 
diffusion mode zd- In the opposite, Q — > limit, they 
both approach their hydrodynamic values: Zh{Q — > 0) = 
DteQ 2 and z D (Q — * 0) = D E Q 2 , with D T e and D E the 
thermal diffusivity and the self diffusion coefficient of the 
Enskog's theory. At intermediates Q values, zt oscillates 
around Zd, with a periodicity dictated by S(Q) and d{Q) 
according to Eq. i|6T|l . 

Sears has calculated the moments of the self part of the 
Van H ove scattering function al large wavevectors ijSearsl 
11972ft . Starting from the general case of a velocity inde- 
pendent central force field, he specialized the result to 
the hard sphere case, evaluating the leading correction 
to the impulse approximation due to final state interac- 
tions, which turns out to be Q : 
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2fe fl m2 A / 0.27 ^ 
m V QIe 



(70) 



The transition from the Fickian to the Gaussian regime 
occurs in this case at QIe w 1- 

Finally, it is worth to recall here one the most sig- 
nificant achievement of the mode coupling theory ap- 
plied to hard sphere fluids, but which has been shown 
to apply t o the wider class of simple fluids. Ernst and 
Dorfman ijErnst and Dorfmanl Il975|) have shown how 
Ea. (|65|l . which retrieves the hydrodynamic expression of 
the sound velocity and attenuation, is actually the lead- 
ing term of an expansion of the kind: 



z h (Q) = a h Q 2 -f3 h Q 5/2 
Z ±{Q) = ±ic Q + a ± Q 2 - 



H. The ionic plasma 



-0(Q 11/4 ) 
(±*-l)/3 ± Q 5 / 2 



0(Q 11/4 ) 
(71) 



The dynamical descriptions given up to this point are 
extensions of models holding for ordinary fluids which, in 
some cases, are modified ad hoc to account for the high 
thermal conductivity of liquid metals. 

A totally different approach is to look at liquid metals, 
from the very start, as a one component plasma, i.e. as an 
assembly of identical, point-like charged particles (ions) 
embedded in a uniform backgr ound (the electrons ) whic h 
neutralizes the total charge l)Baus and HansenL Il980|) . 
The long-ranged, Coulomb interaction active in this case 
give rises to peculiar phenomena such as screening effects 
and plasma oscillations. 

The ionic number density p(r, t) and the current den- 
sity j(r, t) are related through the continuity equation: 



dp(r,t) 
dt 



-V-j(r,t) 



The Poisson equation, for a fluid with Ze charge, reads 



V • E(r,i) = A7rZeSp(r, t) 

in which <5p(r, t) is the deviation of the ionic density from 
its average value p (neutralized by the opposite uniform 
electronic density). Neglecting thermal conductivity ef- 
fect (collisionless regime), from the Newton's law for the 
equation of motion one can find a third equation and close 
the system. In the limit of large wavelengths fluctuations 
(compared to the ionic size) one can write 



m - 



dt 



pZeE{v,t) 



being m the ionic mass. Combining the previous equa- 
tions one easily gets: 



d 2 p(r, t) 
~~dt 2 



which describes ionic plasma oscillation of characteristic 
Q independent frequency 



AirpZ 2 e 2 



(72) 



This simplified picture, therefore, contrasts with the 
experimental evidence of long wavelength excitations 
whose frequency vanishes in the Q — > limit. The com- 
monly used remedy is to " dress" the plasma frequency ac- 
counting for the electron screening effect, i.e. to take into 
account that the background electrons have their own dy- 
namics that can be described introducing the dielectric 
response e(Q). In such a way the coulomb interactions 
and, in turn, the plasma frequency are renormalized lead- 
ing to the expression: 



3_ 

<Q) 



(73) 



In the small wavevector limit the Thomas-Fermi 
expression within the Random Phase Approximation 
(RPA, in with the system is assumed to have a free 
electron gas compressib ility) yelds the useful expression 
l|March and Tosilll99l[) 



<Q) = i 



Qtf 



(74) 



in which Qt f 



h 2 (3ir 2 p e ) 2 / 3 /2m, 



- 6ire pe/Ep, with Ep = 
p e )"' " l'lTn e , p e and TO e the Fermi energy, the 
electronic density and the electronic mass, respectively. 
At small Q values, Eq.dJ, together with Eqs.O and 
(J73J, leads to the so called Bohm-Staver expression, i.e. 
a dispersive excitation with sound velocity 
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(75) 



in which vf the electron velocity at the Fermi level. 

As we will show in the following, the experimental val- 
ues for sound velocities in conductive liquids often con- 
trast with the prediction of Eq.J7SJ, especially at small 
electron density or, equivalently, at large values of the re- 
duced ionic radius r s — (-. — ^-r) 1 / 3 , where the Thomas- 
Fermi approximation is no longer valid. 

The interplay between plasma oscillations and sound 
waves can be better accounted for within a two 
component plasma description (TCP) in which nucle i 
and electrons are treated separately l)Chiharal Il985fk 
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eventually exploiting the framework of Mori-Zwanzig 
ijHansen and Sjogren! ,ll982j) . The details of the TCP are, 
however, beyond the scope of this review, the interested 
reader might want to consult more specialized literature. 



III. EXPERIMENTAL STUDY OF THE MICROSCOPIC 
DYNAMICS 

The main experimental ways to study microscopic dy- 
namics in liquid metals are acoustic spectroscopy and 
inelastic scattering experiments. These latter have to be 
necessarily performed with probes which can penetrate 
enough into the sample to give boundary free informa- 
tion, a requirement which restrict the choice to neutrons 
and X-rays only. Actually, a few attempts have been per- 
formed by means of visible light scatterin g (iDill. 1 1 982t> . for 
insta nce on liquid mercury and gallium l|Dil and Brodvl 
1976), but Brillouin scattering from opaque, liquid mat- 
ter, presents several difficulties, mainly associated with 
the ill definition of the exchanged momentum in absorb- 
ing media. 

The basics of the neutron inte raction with matter 
have been surveyed in many pap e rs llCoplev and Loyesevi 
19751) and bo oks <|Egelstafi Il967t iLovesevl 



Marshall Il97lj) which provide exhaustive survey of this 



issue. Inelastic X-ray scattering is a relatively newer tech- 
nique, and therefore we will detail the basics of an IXS 
experiment recalling from time to time INS features for 
comparison. 



A. The scattering problem 

The measured signal in an inelastic scattering experi- 
ment is determined by the double differential scattering 
cross-section. Within the linear response theory, where it 
is assumed that the coupling between the probe and the 
system is weak, this scattering differential cross-section 
can be written quite generally as the product of three 
terms: i) One term describes the intensity of the probe- 
sample coupling, and it is independent from the energy 
of the incident particle, ii) A second one is a kinematic 
term related to the phase-space volumes of the incident 
and scattered particles, iii) The third term is the space 
and time Fourier transform of the correlation function of 
the observable in the system that couples to the probe 
particle. This last quantity is the one related to the ele- 
mentary excitations characteristic of the system. 



1. The photon-electron interaction Hamiltonian 



(i.e. to zero order in the electron-to-nuclei mass ratio 
m e /m ), and neglecting the magnetic terms (i.e. to the 
zero order in the electron spin) one gets: 



H 



INT = 



2m e c 



2m P c 



rr(l) , rr(2) 
n INT ~r 11 INT 



where the symbol {} denotes the anti-commutator opera- 
tor. The two "electric" terms, Hj NT and Hj 2 2 T contain, 
respectively, two and one field operators A(r). It is clear, 
therefore, that - in a perturbation expansion treatment 
of the interaction Hamiltonian - the term Hj NT will give 
rise to two-photons processes at first order while the term 

(2) 

Hint wm S^ ve r ^ se to one-photon processes at first or- 
der. To have the two-photon processes from the latter 
(the so called p • A contribution), necessary to describe 
the scattering events, one must consider the second or- 
der in the perturbation expansion, which is consequently 
completely negligible in the off-resonance case. In the fol- 
lowing, therefore, we will consider only the first charge 
scattering term. 



2. The X-ray scattering cross-section. 

The double differential cross-section, aE g n , is pro- 
portional to the probability that an incident particle is 
scattered with a given energy and momentum variation 
within an energy range AE and a solid angle AQ. In 
the process, a photon of energy Ei, wavevector ki, and 
polarization e^, is scattered into a final state of energy 
Ef, wavevector ky, and polarization ej, and the electron 
system goes from the initial state |J) to the final state 
\F) (states with energies Ei and Ep, respectively). Ac- 
cording to this definition, the double differential cross 
section can be related to the quantity dF% d ^ f which is the 
probability rate per sample and probe unit that a probe 
particle makes the transition from the initial state to the 
final state: 



d 2 a dFU/1 d 2 n 



dEdVL 



dt j dEdfl 



In this equation j is the incident particle current den- 
sity (j = P v > being p the particle density and v its ve- 
locity) and g E 2 n the density of states of the scattered 



particle. For zero mass particles, the latter two quanti- 
ties can be written as: 



The actual expression for the scattering cross section 
can be derived by a perturbation expansion from the 
probe-system interaction Hamiltonian. In the case of the 
interaction of charges with the electromagnetic field, in 
the weak relativistic limit (i.e. to first order in v 2 /c 2 ), 
neglecting the direct coupling of the field with the nuclei 
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(76) 



(77) 



dEdn 

Therefore, the double differential cross section becomes 
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dEdn 8tt 3 he 2 dt 

The transition of the incident particles between states 
i and / involves, in general, different possible elementary 
excitations in the sample. This implies that, indicating 
with F%,I l jT t ' F ; the scattering probability involving the 

transition in the sample from the state \I) to the final 

dP^f 



state \F), the total probability 



dt 



can be expressed 



dPt. 



dt 



E 



dPi 



dt 



Equation (|7S|) is particularly useful, as the transition 

probability per unit time dP ''' d ^ f ' F can be calculated from 
the perturbation theory. To first order this quantity is 
written as (Fermi golden rule): 



the following considerations: i) One assumes the validity 
of the adiabatic approximation, and this allows to sepa- 
^g-j rate the system quantum state \S) into the product of an 
electronic part, \S e ), which depends only parametrically 
from the nuclear coordinates, and a nuclear part, \S n ): 
\S) = \S e )\S n ). This approximation is particularly good 
for exchanged energies that are small with respect to the 
excitations energies of electrons in bound core states: this 
is indeed the case in basically any atomic species when 
considering values in the range of phonon energies. In 
metals we neglect the small portion of the total electron 
density in proximity of the Fermi level, ii) One limits 
to consider the case in which the electronic part of the 
total wavefunction is not changed by the scattering pro- 
cess, and therefore the difference between the initial state 
| J) = \I e )\I n ) and the final state \F) = \I e )\F n ) is due 
only to excitations associated with atomic density fluc- 
tuations. Using these two hypothesis we then obtain: 
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Inserting the term H^ T into Eq. l(7§|) . using Ea. (|7Sj) . 
and considering the initial and final photon states as 
plane waves one gets: 
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where Q — kj — k/ {E = Ef — Ei) is the momentum 
(energy) transferred from the photons to the system. The 
sum over the initial and final states is the thermodynamic 
average, and Pi corresponds to the equilibrium popula- 
tion of the initial state. 

Apart from the sum over the phase factors of the pho- 
tons scattered from the different particles, whose inter- 
ference give rise to a truly Q dependent scattering sig- 
nal, the energy- and angle- integrated cross section is of 
the order of the square of the classical electron radius 

ni = — ~ — j ■ 



3. The adiabatic approximation and the dynamic structure 
factor 

From Eq. (|80J) . which implicitly contains the correla- 
tion function of the electron density, one arrives at the 
correlation function of the atomic density on the basis of 
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where fj(Q) is the atomic form factor of the j-th atom 
at Rj and the sum is now extended to all the atoms 
(molecules) of the systems. Assuming that all the scat- 
tering units in the system are equal, this expression can 
be further simplified by the factorization of the form fac- 
tor of these scattering units, and by the introduction of 
the dynamic structure factor S(Q, E) defined as: 



S{Q,E)= PiJ{E-E f + Ej)) 



F n J n 



(F v 



3 



By representing the <5-function above as a time integral, 
indicating by (...) the thermal average (o) = J2i ^M-^l-O 
and using the completeness operator \F n }(F n \ = 1 
the dynamic structure factor can be also written in the 
more familiar form: 



S(Q,E) 



2ttN 



3,k 



where N is the number of particles in the system and 
the sum over (j, k) extend over these N particles. The 
double differential cross-section can then finally be re- 
written as: 



d 2 o 
dEdn 



= (- 



-L( ei -e f ) 2 \f(Q)\ Z S(Q,E) (82) 



In the limit Q — > 0, the form factor f(Q) is equal 
to the number of electrons in the scattering atom, i.e. 
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f(Q) = Z. For increasing values of Q, the form factor 
decays almost exponentially with a decay constants de- 
termined by the size of the radial distributions of the 
electrons in the atomic shells of the considered atom. At 
Q-values large with respect to the inverse of these di- 
mensions, therefore, the inelastic X-ray scattering from 
density fluctuations is strongly reduced. 

The cross-section derived so far is valid for a system 
composed of a single atomic species, This derivation, 
however, can be easily generalized to molecular or crys- 
talline systems by substituting the atomic form factor 
with cither the molecular form factor, or the elementary 
cell form factor respectively. The situation becomes more 
involved if the system is multi-component and disordered. 
In this case the factorization of the form factor is still 
possible only assuming some specific distribution among 
the different atoms. In the limiting case that such distri- 
bution is completely random, an incoherent contribution 
appears in the scattering cross-section. 

4. From cross section to count rate. 



p w Z 4 with modifications at energies close to electron 
absorption thresholds. Consequently the effective scat- 
tering volume is very much reduced in materials with a 
high Z (as Z 4 ), while the cross section increases as Z 2 , 
making the study of these materials by all means more 
difficult than for those with low Z. The behavior of the 
optimal signal intensity as a function of atomic number in 
monatomic systems with sample length L = — can be de- 
duced from the data reported in Fig. ^ There, we show 
the quantity ^jf-, with a c — (rgZ) 2 , which gives directly 
a measure of the efficiency of the method at the consid- 
ered photon energy: in this example we took an incident 
photon energy of 22 keV. The quantity ^jf- is by defini- 
tion the ratio between a c and at, where at is the (mea- 
sured) total X-ray cross-section of the considered atom. 
This analysis is useful however only when it is possible 
to study samples of optimal length. In cases where the 
sample size is limited either by its availability or by the 
sample environment (extreme pressure, high/low temper- 
ature, high magnetic field...) it is obvious that one has 
great advantages in studying high Z materials. 



The Z-dependence of the Thomson scattering cross- 
section seems to imply a facilitation in studying systems 
with high Z . In reality, this is no longer true when the 
effect of photoelectric absorption is taken into consider- 
ation. Indeed, neglecting multiple scattering events, the 
signal detected in an IXS experiment from an infinitesi- 
mal slab of thickness 5x orthogonal to ki can be written 
as: 



dN = N {^^)AEAn P dx (83) 

where No is the flux of the incident photons, N is the 
flux of scattered photons in an energy interval AE and 
in a solid angle Ail, p is the number of scattering units 
per unit volume and p is the total absorption coefficient. 
Dealing with a macroscopic sample of length L, in the 
relevant case of a nearly forward scattering geometry, 
Eq.JHSJ becomes 

dN(x) = N e-^^J^-)AEAnpdxe~ fi( - L ~ x) (84) 
which, integrated over the whole sample length, yields: 



N = Noi -dEdn )AEAnpLe ~ (85) 

Let us discuss the L-dependence of this function. It 
is obvious that N attains a maximum (the optimal sam- 
ple length) when L = /U -1 , and that the value of N at 
this maximum point is proportional to pj~ . Considering 
an X-ray energy of approximately 20 keV and Z > 3, 
p is almost completely determined by the photoelectric 
absorption process. This process gives approximately 
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FIG. 1 Ratio between the total number of photons scattered 
by the Thomson process and those lost through all the other 
processes, among which photoelectric absorption, in a sample 
of length /i -1 calculated as a function of the atomic number 
Z for photons of incident energy of ~ 22 keV 

Equation (|85|l accounts for one scattering events only. 
An estimate for the two scattering process intensity can 
be obtained by invoking the forward scattering approxi- 
mation again, indeed, the ratio of the two over the one 
scattering rates 

N (2) / N (l) reads . 

N^ _ nphrUolf(8)S(e)] 2 d6 
NW Z 2 S(0) 

in which h is sample transverse dimension traversed 
by the incident beam. The integral accounts for all the 
possible two scattering process leading to a final forward 
scattering. This expression shows how to suppress multi- 
ple scattering one has to reduce the transverse beam di- 
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mension. A similar estimate for neutron scattering is pre- 
vented by the much more complicated scattering paths, 
since, within a similar Q range, the scattering angle in 
INS is normally much lager. However, it can be noted 
that in the case of neutrons the typical transverse beam 
size is much larger (as 10 -j- 100 mm) than the IXS ones 
(« 100 /im), thus resulting in a more important contri- 
bution requiring accurate corrections. 



5. Kinematics of the scattering processes. 

Another important difference between X-rays and neu- 
trons scattering lies in the kinematics of the scattering 
processes. The momentum and energy conservation laws 
impose that: 
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Q = k, k/ 

E = Ef — Ei 



Q 



kt — 2kikfcos0 



where is the scattering angle between the incident 
and scattered particles. The relation between momentum 
and energy in the case of photons is given by: 



FIG. 2 Kinematic region accessible to IXS in reduced E/Ei 
and Q/ki units (left panel). The right panel shows the real- 
istic cases of molten lithium (open circles) and sodium (full 
circles) for an incident energy of 25 keV. 



and, therefore: 



E(k) = hck 
and therefore one obtains: 



h , 



I*) - 2 (l- f] corf 



Considering that the energy losses or gains associated 
to phonon-like excitations are always much smaller than 
the energy of the incident photon (E « Ei) this relation 
reduces to: 



2sin— 
2 



(E « Ei) 



This last relation shows that, in the limit of small en- 
ergy transfers, the ratio between the exchanged momen- 
tum and the incident photon momentum is completely 
determined by the scattering angle. Therefore, in inelas- 
tic X-ray scattering, there are basically no limitations in 
the energy transfer at a given momentum transfer for 
phonon-like excitations. 



6. X-rays vs. Neutrons. 

At variance with the previous equation, if the probe 
particles have mass m p 



E(k) = 



2 1.2 



K z k 
2m r 
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In this case of thermal neutron scattering, the approx- 
imation E << Ei no longer holds, and the kinematics 
of the scattering experiments is determined by Eq. 187|l . 
As an example, in Fig. |3| we report the accessible kine- 
matics regions in the |- vs p plane for two different 
incident energies, indicating paths at constant scattering 
angles. In the same figure, similarly to Fig. |3 we also re- 
port the approximate dispersion curves for liquid lithium 
and sodium. In the best situation, i.e. in forward scat- 
tering where the accessible region is as wide as possible, 
the limiting curve is linear around Q = 0(E = 0), and 

its tangent is M: — 2^-. Recalling that Ei = -^""S onc 
gets E = vn?iQ, with vn the velocity of the incoming 
neutron. As the dispersion relation for acoustic phonon 
is linear, E = v s HQ, with v s the velocity of sound, it 
is clear that whenever v s is larger than vn the excita- 
tions peaks lie outside the accessible region and, there- 
fore, when v s > vn the neutron technique cannot be 
applied to study the acoustic branch. This limitation 
does not apply to the case of X-rays, as, according to Eq. 
(|86|l . there are basically no limits to the energy region 
accessible at a given scattering angle. 

As discussed before, the presence of relevant absorp- 
tion phenomena is the main effect that determine the 
scattering volume in an IXS experiment. This implies, 
therefore that the probability that a photon is scattered 
from the sample is small, and this strongly suppresses 
the multiple scattering processes. In IXS experiments, 
indeed, the multiple scattering can be disregarded, thus 
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while in the X-rays case kt ss fej , and this factor is 
very close to 1, in the neutron case this term give 
rise to a Q dependence of the scattered intensity. 

3. No polarization terms are present in the cross sec- 
tion for neutron, while in the case of X-rays the 
term tells us that the Thomson scattering arises 
from a scalar interaction and therefore the polar- 
izations of the incident and scattered photons must 
be parallel. 
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FIG. 3 Kinematic region accessible to neutron scattering ex- 
periments for incident energies Ei = 50 meV (left panel) and 
Ei — 500 meV (right panel), reported for different scattering 
angles. Open and full circles are the (approximated) sound 
dispersion of molten lithium and sodium, respectively. 

avoiding the use correction procedures. This is there- 
fore an important advantage with respect to the neutrons 
case, where, on the contrary, the sample length is deter- 
mined by the scattering (rather than absorption) length. 
In Fig. 0|are reported for comparison the (coherent) scat- 
tering lengths of the elements for the X-rays and neutrons 
cases. 

Finally, it is worth to compare the double differen- 
tial scattering cross section for X-rays obtained before 
(Eq. (|82|) ) with the similar quantity derived for neutron 
scattering (in the hypothesis of fully coherent scatter- 
ing). The latter quantity is not derived here for brevity 
and can be found in many te xtbooks ijEgelstafl Il967t 
lLovesevt Il987t iMarshallL Il97l|) . The two cross-sections 
read as: 

d 2 = f r 2 ^ (e^) 2 \f(Q)\ 2 S(Q,E) xray 
dEdn \ b 2 % S{Q,E) neutron 

(88) 

Beside to the proportionality of both the cross sections 
to the dynamics structure factor, it is worth to underline 
that: 

1. The two cross section are proportional to a char- 
acteristic scattering length squared (ro in the case 
of X-ray and b in the case of neutrons) that are 
comparable in magnitude (see Fig. ). The further 
factor \ f(Q)\ 2 (oc Z 2 for small Q's) in the case of X- 
ray does not bring an increase of the actual signal 
in the experiment because, as discussed before, in- 
creasing Z also limit the optimal scattering volume 
due to the increase of the photoelectric absorption. 

2. In both cases the phase space of the incident and 
final plane waves gives rise to the factor , however 



4. Finally, the X-ray scattering cross section contains 
the form factor f(Q), i.e. the Fourier transform 
of the charge density spatial distribution. As the 
charge density is localized around the nuclei in a 
space region of typical dimension of few tenth to 
few hundredth of nm, the function f(Q) decreases 
appreciably on a Q range of several inverse nm -1 , 
thus it does not depresses too much the scatter- 
ing cross section in the mesoscopic region of inter- 
est. In the case of neutrons this form factor is not 
present (actually it is equal to 1) as the neutrons in- 
teract with the nuclear matter, localized in typical 
dimension of 10~ 6 nm. The neutron form factor is 
therefore constant in the whole accessible Q region. 

From the discussion made so far, it should be now quite 
easy to understand how important is the development of 
the X-ray method, which can access, in principle, an ex- 
tremely large region of the E — Q plane. Particularly 
important is the small Q region, where the acoustic exci- 
tations have energies which are not of easy access to the 
neutron spectroscopies. 
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FIG. 4 The coherent scattering cross section of the elements 
for X-ray (open circles) and for neutrons (full circles) reported 
as a function of the atomic number Z. 



B. From the experimental data to the dynamical quantities 

In order to extract quantitative information from the 
experimental intensity, i.e. to perform measurements of 
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S q (Q,uj) on an absolute scale, the most direct way is to 
use a reference scatterer and this is customarily done in 
neutron experiments. In IXS, for instance, such a proce- 
dure can be quite difficult because of the Q-dependence of 
the form factor and of the analyzers efficiencies. For these 
reasons indirect method are always preferred. One possi- 
bility is to explo i t the lo west order sum rules of S q (Q, u>) 
ijScopigno et all |20 00b): in particular for the first two 
frequency moments one has the 

(cu a ) Sq = J S q {Q,u)dw = S{Q), 
(uj 1 )s q = J u;S q (Q,u})duj = hQ 2 /2m. 

where the second equality follows from Eq. (JSJ applied 
for n — 1. The measured raw intensity is related to the 
dynamic structure factor through 



I(Q, u) = A(Q) J dtu'SgiQ, lu')R(lu - J) (89) 

where R(lo) is the experimental resolution function and 
A(Q) is a factor taking into account the scattering geome- 
tries, the experimental setup and the atomic form factor. 
The first moments of the experimental data, (u>°)i and 
(w 1 )/, and those of the resolution function, (cu°)r and 
are related to (cj°)s and (^ 1 )s by: 

(^j = A(Q)(u> ) Sq (u,°) R , 

(a, 1 ), = A(Q){(u, a ) St {u, l )R + (w l )s,{u ) R ). 

From the previous equation one derives that 



S Q (Q) = 2^((^)i/(^°)i - {^r/^r)- 1 . (90) 

This procedure, therefore, can been adopted to estab- 
lish an absolute scale for S q (Q,tu) using the experimen- 
tally determined I(Q,lu) and R(uj). 

C. Handling liquid metals 

Working with liquid metals poses several practical 
problems. In particular alkali metals are highly reactive 
and need to be kept under a protective atmosphere. A 
relatively small impurity (less than 100 ppm) in Ar or 
nitrogen will cause a film to form on the surface of the 
liquid metal. 

In addition, glass is often limited in its use as a con- 
tainer for most liquid metals. Liquid metals are often 
strongly reducing. Glass, composed chiefly of silicon 
dioxide, is penetrated by the metal atoms, which can 
reduce the silicon by forming a metal oxide. As a result, 



the glass becomes discolored and brittle. For these rea- 
sons Pyrex can not be used above 600 K and pure quartz 
above 900 K. 

Preferred materials for working with liquid metals are 
the refractory metals. This refers to the titanium group 
(Ti, Zr) as well as the vanadium and chromium groups. 
These transition metals are much less likely to undergo 
reduction and be solvated by liquid metals. The disad- 
vantage in addition to cost is that the refractory metals 
have certain properties that make fabrication difficult. 
A trade-off is to use iron plated with chromium. Other 
less reactive transition metals, such as the noble metals, 
are soluble in many liquid metal solutions. Austcnitic 
stainless steel can be suitable up to 1000 K 

Liquid metals also share a common chemistry. The 
increasing electropositivity of the metals composing the 
liquid metal solution will determine the liquid's reactiv- 
ity. Mercury, which is not very electropositive, is stable 
in air. Alkali metals, which are the most electropositive 
group of elements, are air-sensitive. Li reacts slowly with 
air, yet dissolves and reacts quite readily with nitrogen. 
The other alkalis are insensitive to nitrogen but react 
with other gasses. All the alkali metals violently react 
with water on this basis: Cs > Rb > K > Na >> Li 
l|Ohselll985|) . The hydrogen generated in the water-alkali 
reaction can, in turn, react explosively with oxygen. 

Reactivities reflect those of the solid material. The 
more reactive liquid metals, usually alkali ones, often de- 
velop a film if exposed to air. The high surface tensions 
of many liquid metals enable this film to remain in place. 
However, rates are greatly increased if turbulence is in- 
troduced because the protective coating is often dissolved 
into the liquid metal. Products of reaction, such as hy- 
drides and oxides, are often redissolved into the liquid 
metal solution. 

By virtue of the previously mentioned difficulties, ex- 
periments with alkali metals are often very hard to per- 
form. In the case of Inelastic X ray Scattering, in partic- 
ular, such difficulties are enhanced by sample dimension 
requirement. As we have seen in sec. I111.A.4I indeed, in 
an optimal IXS experiment the sample length has to be 
comparable with the absorption length. With a few ex- 
ceptions (Li and Na), this typically means sub- millimeter 
sample thickness. Common choices are, therefore, sam- 
ple cells made of compatible metals provided with sealed 
sapphire or diamond windows. 

Large efforts have been done recently to overcome the 
difficulty of performing X-ray experiments on liquid met- 
als, the most remarkable example being the so called 
Tamura-type cells made with a single c rystal sapphire 
with Be windows pressurized under He ijTamura et all 
1999), performing up to 1900 K and 2 Kbars. More re- 
cently, a new sample environmen t especially tailored fo r 
alkali metals has been proposed ijMatsuda et all [2004) . 
In this case the cell is entirely made of molybdenum, with 
the windows made by single- crystal disks of controlled 
orientation electrolytically thinned at w 40 /im. 

A totally different approach is the one of contact-less 
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FIG. 5 Coupling IXS with levita tion techniques: C onstant 
spectra of liquid alumina. From (ISinn et ttZ.Ll2003) 



techniques. In this case the sample is levitated either 
electrostatically or by means of a controlled gas jet. The 
main difficulty in this technique is related to the sample 
stability in the x-ray beam, but recent impressive ad- 
vancement have been done in this field. A new beamline 
for electrostatic levitation (BESL) has been developed at 
the Advanced Photon Source (APS), and the relevance 
of icosahedral ordering in the superc ooling capabi l ities o f 
liquid metals has been investigated ijKelton et adl2003|) . 
Another example is a recent X-ray scattering experiment 
performed on liquid AI2O3 in which alumina droplet of 3- 
4 mm diameter have been levitated b y gas jet flow on t he 
inelastic scattering beamline 3ID-C ( Sin n et "oil I2003) . 



IV. EXPERIMENTAL RESULTS 

In this section we review, to the best of our knowledge, 
the experimental results reported so far, ordered accord- 
ing to the sample group in the periodic table. No results 
are available so far for elements belonging to group II. 



Alkali metals do not occur freely in nature, they are 
very reactive and can explode if exposed to water. These 
metals have only one electron in their outer shell and, as 
with all metals, they are malleable, ductile, and good 
conductors of heat and electricity. Alkali metals are 
softer than most other metals. Among the metallic el- 
ements they share the simplest pairwise interaction po- 
tential, which is also the closest to the Lennard Jones 
one. As a consequence, their structural properties are 
also particulary simple, with a structure factor resem- 
bling the one of hard spheres. Also the dynamics, there- 
fore, is expected to mimic the theoretical and numerical 
results achieved for Lennard Jones and hard sphere sys- 
tems. 



1. Lithium 

Liquid lithium is probably the system which better re- 
veals the complementarity of neutrons and X-rays as far 
as inelastic scattering is concerned. Due to the high ab- 
sorption cross section of the 6 Li isotope (a a = 940 b) 
neutron scattering experiments must necessarily be per- 
formed on 7 Li enriched samples, which is the dominant 
specie in the natural abundance. The high sound veloc- 
ity (ct ~ 4500 m/s), and the almost equivalent neutron 
scattering cross sections (<jj = 0.68 b; a c — 0.62 b for 
7 Li), pose severe limitations to the use of INS aiming 
at the determination of collective properties, while this 
technique turns out to be extremely useful for the in- 
vestigation of the single particle motion. The first INS 
studies on this system c an be traced back to the work 
of De Jong and Verkerk llDe Jond . Il993t iDe Jong et all 
1993: IVerkerk et all 1 1 992(1 . who showed the presence of 
collective modes with a series of accurate experiments. 
Though they had to face the above mentioned draw- 
backs, indeed, they were able to point out some sig- 
nificant issues: by modelling the coherent contribution 
with the extended hydrodynamic model (see Eq. (|63|l of 
section III.Gj l they measured the dispersion curve above 
Q m /2 (see Fig. [7J and they reported deviations from 
the Landau-Plazek ratio, which is expected to hold in 
the hydrodynamic regime (see Eq. H18|) of section lll.Cfl . 
On the other side they shed light on the single particle 
motion, accurately determining the incoherent contribu- 
tion to the dynamic structure factor within the frame- 
work of section III.GI (Eqs 12813 II)'). They corroborated 
the M ode Coupling predictions (|de Schepper and Ernstl 
119791) . extracting values of the diffusion coefficient and 
determining its temperature dependence. 

An exhaustive characterization of the coherent dynam- 
ics was provided by the advent of Inelastic X ray Scat- 
tering developed in the early nineties, and liquid lithium 
has been the benchmark of such development. Being the 
lightest of the liquid metals, indeed, lithium played a 
privileged role in IXS, for the favorable signal to noise 
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to exploit IXS to gather insight into the microscopic dy- 
namics of disordered systems. 



FIG. 6 Pioneering (1991) low resolution IXS determination of 
the dy namic structure factor in liquid lithium with INELAX. 
From l|Burkel and Sinnl Il994ft . 
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FIG. 7 Dispersion curve of liquid lithium achieved with IXS 
at INELAX. Theoretical pred ictions and INS results at higher 
Q's are also reported. From iBurkel and Sinnl 11994) . 



ratio and for the high sound velocity which allowed to 
resolve the inelastic spectral component minimizing the 
initial difficulty of achieving energy resolutions compara- 
ble to neutrons. 

Since the pioneering work of Burkcl llBurkeli 11991ft 
with the INELAX instrument (see fig. |S] and [JJ , a de- 
cisive step forward achieved with the advent of the third 
generation sources which, combined to a brilliant tech- 
nique for manufacturing silicon crystal analyzers, allowed 
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FIG. 8 First IXS measurements on li quid lithium per formed 
on a third generation source (ESRF) (ISinn e t at, 1997[). En- 
ergy resolution is here 5E = 11 meV. The continuous line is 
the best fit accord ing to the extend ed hydrodynamic model 
of Eq. E3. From JSinn et qZ.Ul997ft 

The first remarkable result on a third ge neration fa- 
cility (ESRF) was provided by Sinn et al. ijSinn et all 
Il997ft who, measuring energy spectra at fixed wavevec- 
tors, reported clear evidence of collective modes, be- 
ing able to give significant hints for the choice of the 
most appropriate pseudopot ential to describe liqu id met- 
als in numerical simulations l)Canales et aZJ . Il994ft . In the 
same work, following the extended hydrodynamic model 
outlined in section J1LGI ijde Schepper and Cohenl Il980t 
Ide Schepper et all Il983ft . it was also reported evidence 
for positive dispersion, i.e. for a sound velocity value ex- 
ceeding the hydrodynamic one. This phenomenon was 
ascribed to a transition from a liquid to a solid-like re- 
sponse. 

Following the development of the IXS technique, new 
experiments have been more recently performed on liquid 
Li in the extended region 1.4 to 110 nm" 1 , corresponding 



28 



to Q/Qm ~ 5 • 10 2 -j- 5, which are reported in Fig. 
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FIG. 9 IXS measur ement of liquid lithium in a wide energy 
momentum- region IjScopjgno et all fcOOOaf) . The transition 
from hydrodynamic to gaussian-like response (continuous line 
in the right panel) can be clearly noticed. Energy resolutions 
are here 1.5, 3.0 and 7.0 meV, increasing with the exchanged 
momentum 

In Fig. El is reported the dispersion relation deter- 
mined in the same energy-wavevector region, and the 
transition between the two distinct dynamical regime is 
here evidenced by the sound velocity behavior. Beyond 
the first quasi-hydrodynamic region (an initial nearly lin- 
ear dispersion), structural effects take place suppressing 
the sound propagation around Q m /2 due to strong neg- 
ative interference. With increasing Q values, the points 
in Fig. El show a second pseudo-BZ, followed by a series 
of oscillations that damp out with increasing Q - here, 
loi{Q) is approaching the single particle behavior. These 
oscillations are in anti-phase with those of S(Q) and are 
therefore associated with the local order in the liquid. 

While at low Q's the dynamic structure factor is qual- 
itatively described by an extended hydrodynamic treat- 
ment (Eq. l|63l of section Hi. G|) . at wavevectors distinctly 
larger than Q m the single particle response is attained 
through the mechanism described in section III.DI ulti- 
mately leading to the expressions well accounted by a 
combination of two Eq. I|24|l . accounting for each of the 
two isotopes 6 Li and 7 Li. 

Thanks to the improvement in the energy resolution, 
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FIG. 10 Sound velocity as deduced by the maxima of the 
current correlation spectra, from the best fit with quantum 
corrected and resolution convoluted models. 
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FIG. 11 Memory function at work: refined lineshape analysis 
of high resolution IXS spectra. Both thermal and viscous 
channel are taken into account, mimicking this latter with 
one (Eq. 156L dotted line) or two (Eq. H59II . continuous line) 
exponential processes. 



which is nowadays comparable to the one of INS spec- 
trometers in the same energy- wavevectors domain (1.5 
meV at present), an approach based on the generalized 
hydrodynamics has been developed, which allowed to 
point out the presence and the role of relaxation pro- 
cesses driving the col lective dynamic at the m icroscopic 
probed wavelenghts ijScopigno et all 2000b c). Within 
a memory function framework IjMorit Il965j) , it has been 
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ascertained the presence of two distinct viscous relax- 
ation channel (see fig. beyond the thermal relaxation, 
clarifying the origin and the nature of sound dispersion 
and attenuation properties in simple fluids. Of the two 
processes, active over well separated timescales, one is 
related to the well known transition between a low fre- 
quency, liquid-like response to an high frequency, solid- 
like response. The second mechanism is instead a general 
relaxation process peculiar of the vibrational dynamics 
which is present regardless the thermodynamic state of 
the system. In this context, the positive dispersion has 
been shown to be strongly related to this latter process, 
being the solid like response already attained over the 
wavevectors range probed in IXS (or INS) experiments. 
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FIG. 12 Constant energy slices of the dynamic structure fac- 
tor determined by IXS. Umklapp modes are visible on the 
sides of the main structure factor peak 

An alternative route to the investigation of collec- 
tive dynamics, which is dual to the one followed in the 
above mentioned experiments and which easily achiev- 
able through IXS, is the determination of the dynamic 
structure factor performing Q-scan for fixed values of 
the energy transfer, reported in fi g. IT21 for the case 
of lithium ijScopigno et all l2001b(l . In this way, one 



is able to have a direct sight over the so called umk- 
lapp modes, i.e. excitations characterized by wavevectors 
which differ by multiple of the reciprocal lattice spacing, 
which have been early reported by means of INS in liquid 
lead llCocking and Egelstafll Il965al iDorner et all Il965l 
iRandolph and Singwilll96rfi . 



2. Sodium 

Pioneering experimental determinations of the scatter- 
ing law in liquid sodium can be traced back to the time 
of the IAEA symposiu m held in Chalk River ijCockingl 
Il963t iRandolphl Il964j) . In this system the ratio between 
the incoherent to coherent cross section is very close to 
one (see table , therefore the separation between the 
two contributions is of crucial importance. Soon after 
Randolph's experiment, his data were an alyzed in terms 
of me an square displacement of an atom ijDesai and Yipl 
1969). This framework (described in section lTl.E.lf) poses 
on the Gaussian assumption for the incoherent cross sec- 
tion, while the coherent contribution was e valuated ac- 
cordi ng to the effective mass approximation l)De Gennesl 
1959): 
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The mean square displacement was then determined 
describing the atomic motion in terms of independent 
harmonic oscillators of frequency luq and lifetime to, 
which, in turns, are related to the spectral density of 
the velocity autocorrelation function }(oj): 
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The basic ingredients of this approach are, therefore, 
the knowledge of the static structure factor, of the macro- 
scopic diffusion coefficient and of the mean squared force 
((VC/) 2 ). The results of this description are tested 
against the experimental data in fig. 1131 From the same 
figure, it emerges how the single particle regime is already 
attained at the lowest reported Q value, i.e. Q = 12 
nm -1 , while the low frequency discrepancy has been ten- 
tatively ascribed to finite instrumental resolution and to 
multiple scattering effect. 

Twenty years la te r, ne w INS data were reported 
llMorkel and Glaserl H 986t HHHi 

rstrom and DahlboreL 
Il984f) . addressing in more detail the incoherent scatter- 
ing contribution and showing how the diffusion process 
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FIG. 13 Randolph's measurement iRandolphl Il964fl of the 

dynamic structure factor of liquid sodium in reduced units: 
a = h 2 Q 2 /2mT and /3 = hu/T for four different val- 
ues of momentum t ransfer. Lineshape analysis according 
to different models iDesai and Yid. Il969f) Continuous line: 
EMA+calculated MSD. Short dashed line: EMA+computer 
simulation computation of the MSD. Long dashed line: hy- 
drodynamic prediction. Dotted-dashed line: free streaming 
limit 
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FIG. 14 Reduced quasielastic linewidth LO1/2/DQ 2 in liquid 
sodiu m at three different temperatures iMorkel and Glaserl 
Il98€f) . The dotted line is the Fickian limit while dash dotted 
line is the perfect gas behavi or (oc 1/Q). The dashed line 
is the hard sphere prediction ( Cohen et all I"l987). while the 
continu ous line is the result obtai ned within mode coupling 
theory iGotze and ZippeliusLll976l) 



is actually more complex. Morkel and Glaser, following 
for the coheren t contribution the Lovesey's prescription 
l|Lovesevl 1 1 9 7 ll) . and ad opting for the incoherent part 
the Nelkin-Gata k mode l ijNelkin and Ghatakl. \1964) de- 
scribed in section ITl.E. 41 extracted the reduced halfwidth 
UJ1/2 of the incoherent contribution finding a crossover be- 
tween the hydrodynamic (lorenzian) and single particle 
(gaussian) regimes. In fig. [Hthe linewidth uj 1/2 /DQ 2 
is reported, and it clearly emerges how the single par- 
ticle limit (wi/2/DQ 2 oc 1/Q) is not yet attained even 
at Q « 40 nm _1 , which contrast the earlier assumptions 
of Desai and Yip. After a low Q diffusion retardation 
the mobility increases in the transition region and finally 
tends to the free gas limit. The whole Q dependency 
is well descr ibed within the Enskog's hard sphere gas 
l|Searsl IT972(1 . in terms of the expression lf7U|) . An alter- 
native description of the single particle dynamics can be 
recovered within the memor y function approach, thoug h 
it fails in the high Q region ijGotze and Zippeliuslll97q) . 

The first IXS determination of the collective dynam- 
ics in liquid sodium is due to Pilgrim and collaborators 
(Pilg rim et aZJ .1 1999) . In this work, the coherent dynamic 
structure factor was measured at several temperatures, 



and analyzed according to the extended hy drodynamic 
mode l previously applied in liquid lithium flSinn et all 
Il997j) . The take-home message is the presence of a pos- 
itive dispersion effect which does not show significant 
temperature dependence. This result seems to rule out 
an interpretation of the positive dispersion in terms of 
an activ ated process, as is the case in hydrogen bonding 
systems (|Monaco et all \l99$) . 

IXS experiments on liquid lithi um were then repeate d 
with increased energy resolution l|Scopigno et all 12002a). 
and analyzed within the same two vi scous relaxation 
proces ses proposed for liquid lithium ijScopigno et all 
l2000cl) . The same data have been also interpreted within 
the fra mework of the scale iny ariance of relaxation pro- 
cesses l)Yulmetvev et all 20031). a theory origin ally de- 
veloped for liquid cesium ijYulmetvev et ad feOQl). which 
has been recently shown to be equivalent to the memory 
function approach in the sense that one solves the chain 
of equations 1351 with some ad hoc closure relation. 
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FIG. 15 Left panel: IXS determination of the S(Q,ui) in liq- 
uid sodium for selected temperatures. Right panel: dispersion 
curves at different temperature. 



3. Potassium 

The first experimental data on liquid potassium ap- 
pea red surprisingly late relatively t o the other liquid met- 
als ijNovikov et all Il996afbl Il997|) . Moreover, the kine- 
matic region Q — E spanned in this experiment was quite 
narrow (10 < Q < 13 imr 1 ) and only partial information 
on the microscopic dynamics could be obtained. 

Very recently, two sets of INS experiments have been 
reported in molten K just a bove the melting tem pera- 
ture, one at the ISIS s ource llCabrillo et all 1200^ . and 
the other at the ILL ijBove et all l2003fl . In the first 
case, two time of flight spectrometers were utilized (IRIS 
and MARI) aiming at a combined study of the dynamic 
structure factor with different energy resolutions for the 
narrow quasielastic and the broader inelastic component. 
The experiment of Bove et al. has been instead per- 
formed on the triple axis spectrometer INI optimized 
to access a broader kinematic region, as shown in fig llGI 
where a detail of the energy-momentum region accessed 
in the two experiments is shown. 

The data taken on IRIS allowed for an accurate, high 
resolution, determination of the diffusive processes un- 
derlying the incoherent dynamics. The results support 
the hydrodynamics predictions corrected by the mode 
coupling terms (see eqs. (|28I31|) of section Hl.E.3|) . as re- 
ported in fig. El Beyond the diffusive mode, Cabrillo 
et al. identify a second contribution, coherent in nature 
according to the authors, which is weaker, broader and 
almost Q independent up to Q Q m , while becomes 
narrower above Q m . The Q dependence of this coher- 
ent mode is rationalized in terms of extended heat mode 
(Eq. (|fci7|> of section lTl.Gfl but its ultimate origin is rather 
ambiguous, especially in view of the INS measurements 
taken at INI and MARI. 
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FIG. 16 Sketch of the kinematic regions accessed in the ex- 
periments of Bove et al. (dashed line) and Cabrillo et al. 
(continuous lines). The linear sound dispersion is also re- 
ported. 
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FIG. 17 Quasielastic linewidth according to recent INS mea- 
surements. Full and open circles are the coherent and in- 
coherent contribution s, respectively, determined with TOF 
ijCabrillo et all l2002f ) . Open tr iangles are the in coherent 
linewidth as measured with TAS iBove et all 120031) . the ob- 
served 2 meV coherent contribution is also indicated. The 
lower dotted line indicates the fickian approximation. 



The results of these two experime nts are repo r ted in 
fig. IT51 for s i milar fixed Q values l)Bove et all 120031 
ICabrillo et all l2002h . As can be easily noticed the pos- 
sibility (offered by INI) of extending the INS measure- 
ments at low Q is paid in terms of resolution. In both 
cases, however, evidence for inelastic coherent scatter- 
ing is reported, though the incoherent scattering largely 
dominates in the region where collective modes are more 
visible. The two sets of data have been analyzed ac- 
cording to different approaches by the respective authors. 
Cabrillo et al utilized a memory function approach trun- 
cating the continued fraction at n = 2, motivating this 
assumption as necessary to account for the nearly Q in- 
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dependence of the coherent quasielastic contribution re- 
ported in fig. El Odd enough, as evinced from Fig. 
1181 (left panel), neither the inverse relaxation time, nor 
the raw quasielastic width that they extract with this 
model favourably compare with the coherent linewidth 
reported in fig. El Bove et al., on the other side, utilize 
a Damped Harmonic Oscillator for the purely inelastic 
term and two lorentian for the quasielastic coherent and 
incoherent contributions, respectively. The results for 
the incoherent part, achieved within the jump diffusion 
model described in section lTl.E.2l are consistent with the 
high resolution measurements (IRIS) of Cabrillo et al. 
(sec fia ll7ll . On the other side, the coherent contribu- 
tion turns out to be much broader (FWHM« 4 meV, 
i.e. a relaxation time r ~ 0.32 ps), in contrast with 
the Cabrillo's data reported in figEI(full circles), but in 
qualitative agreement with the quasielastic linewidth and 
the relaxation time of the same author's measurements 
reported in figfTFI 
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FIG. 18 Left panel: TOF determination o f the DSF in liquid 
potassium (circles) iCabrillo et all [2002f) . The dash-dotted 
line depicts the coherent contribution. Right panel: TAS 
meas urements in a sim ilar momentum transfer region (cir- 
cles) llBove et all l2003f) . The continuous line is the inelastic 
contribution to the collective dynamics 

Both the experiments extract the dispersion curves, 
in one case following the exact hydrodynamic prescrip- 
tion as the maxima o f the current correlat ion function 
llCabrillo et all l2002|) and, in the other ijBove et all 
2003) as the DHO frequency, which coincides with the 
current correlation maximum if the presence of the 
quasielastic coherent term is neglected. The two inde- 
pendent determinations are indeed in good agreement, 
except at large wavevectors where the data of Bove et 
al. are systematically higher though with some scatter- 
ing. The sound velocity values exhibit the usual excess in 
respect to the hydrodynamic value. This high frequency 
sound is ascribed by both the studies as a reminiscence of 
solid like behavior, i.e. as the upper edge of a transition 
occurring from the low Q, hydrodynamic domain to the 
high frequency regime. This claim stems on the basis of 
the similarity of the sound velocity value of molten potas- 



sium with the value for the crystalline acoustic phonons 
along the [1 0] direction. 
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FIG. 19 High reso lution IXS measurem ents in liquid potas- 
sium (open circles) ((Monaco et fflZ.l l2004") . The continuous line 
is the lineshape description according to a multiple relaxation 
memory function model (see text). 



A r ecent IXS experiment on molten K l)Monaco et all 
12004ft contributed to shed some light on some aspects 
of the collective dynamics, giving a coherent picture in 
terms of relaxation processes which is common to several 
other simple fluids and, more generally, to glass forming 
materials and molecular liquids (see fig ll9fl . 

First, it has been shown how the coherent dynam- 
ics is driven by thermal and viscous processes. These 
latter, which are dominant, proceeds over two different 
timescales. Consequently, the FWHM of the quasielas- 
tic (coherent) contribution is per se not directly asso- 
ciated to any relevant timescale. The thermal process, 
indeed, is characterized by a timescale largely exceeding 
the Brillouin frequency, while both the viscous processes 
controls the quasielastic width. The corresponding relax- 
ation times can be instead determined within the memory 
function formalism of Eq. I|59(l , obtaining the results re- 
ported in fig. |2SI In the same plot, the relaxation time 
obtained by Cabrillo et al. (consistent with the deter- 
mination of Bove et al.) is also reported. As one might 
expect, this value is somehow averaging between the two 
distinct viscous processes. More specifically, at low Q 
the one time approximation of Cabrillo et al. seem to 
mimic the slower process, while at higher Q is describing 
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FIG. 20 Viscous relaxation times as measured by means 
of IXS iMonaco et all l2004h . Circles: structural relax- 
ation time. Triangles: microscopic relaxation time. The 
relaxa tion time obta ined by means of INS is also reported 
ijCabrillo et all [2002ft . showing how it averages between the 
two mechanisms reported by IXS. 
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FIG. 22 Sound velocities as determined by IXS: apparent 
(circles, from the maxima of Cl(Q,u>), isothermal (dotted 
line) and infinite frequency limit (continuous line) determined 
from Eas. l|46|l and ilo'il . respectively. The unrelaxed sound 
velocity values, due to the structural relaxation only (down 
triangle) and to the whole relaxation process (uptriangle) are 
also reported as estimated by the fitting. 



4. Rubidium 
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FIG. 21 Dispersion curves (maxima of the cur rent correla- 
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the faster process. This hypothesis is consistent with the 
observation, reported in other alkali metals, of a decreas- 
ing weight of th e slow relaxation proce ss on increasing 
the wavevector l|Scopigno et al 1 12000b). As far as the 
sound propagation properties are concerned, the IXS ex- 
periment analyzed in terms of generalized hydrodynam- 
ics suggest a minor role of the structural process, which 
accounts for approximatively 10% of the whole positive 
dispersion effect, which is dominated by the faster pro- 
cess (see fig. |22J). This observation, already reported for 
many other simple liquids, poses against the commonly 
invoked explanation of the positive dispersion in terms 
of transition from liquid to solid like regime (structural 
relaxation) . 



Liquid Rubidium has been the first of the alkali met- 
als to be ad dressed by a very famous neutron scattering 
experiment l|Conlev and RoweL Il974l) . immedi ately sub- 
stant iated by molecular dynamic simulations ijRahmanl 
11974 . The reason for such interest lies in the possibil- 
ity of extracting information on the collective dynam- 
ics, given the almost negligible incoherent cross section 
(<Tj/(T c « 10 -4 ) and the relatively low sound velocity 
value. The result of this experiment contributed to open 
the route to the understanding of the collective dynamics 
in simple liquids, showing that the presence of an high 
frequency inelastic mode is an intrinsic property of the 
alkali metals not related to quantum properties or critical 
thermal population effect as suggested by earlier works 
on liquid hydrogen ijCarneiro et all Il973t> . Though the 
experiment was affected by an elaborated multiple scat- 
tering subtraction (due to the lack of an absolute normal- 
ization), which lead to a possibly unreliable quasielas- 
tic spectral component, some other chords of interest 
where hit. At varian ce with earlier results on liquid lead 
Iporner et ad Ii967). no evidence of secondary modes of 
transverse nature was reported in this system. Finally, a 
mild positive dispersion effect was observed (though Cop- 
ley and Rowe looked at the maximum of S(Q,lo) rather 
than to the maximum of J(Q,u>)) which was tentatively 
ascribed to the distinctio n of zero sound and first sound 
as discussed by Egelstaff ijEgelstaffl Il967|) . 

More recently, an inelastic s cattering experimen t was 
performed with cold neutrons l)Chieux et all Il996|) aim- 
ing at the determination of the dynamic scattering law in 
an extended temperature region beyond the one explored 
by Copley and Rowe. The kinematic accessed region is 
above Q — 9 nm -1 , and therefore the observed excita- 
tions lies beyond the linear dispersion region. This work 
has the merit to stress the importance of the choice of the 
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FIG. 23 The experimental determination of constant Q slices 
of the DSF by means of INS scattering in liquid Rubidium. A 
new er a for the study of collecti ve properties in simple liquids. 
From (ICoplev and Rowa, Il974f) 



appropriate dynamical variable and of the fitting model 
to determine the dispersion curve. 

The old data of Copley and Rowe have been more re- 
cently reanalyzed in terms of gen eralized hydrodynamics 
ijMorkel and Bodensteinerl ^990) , comparing the results 
to the ones obtained in molten Cesium, which are dis- 
cussed in section ITV.A.5I 

To our knowledge no IXS measurements have been re- 
ported on liquid Rubidium. The main difficulties for such 
an experiment would be the very small absorption length 
(about 200 /im) and the quite low sound velocity value 
(« 1400 m/s) which would confine the clastic modes on 
the tail of the resolution function. 



5. Cesium 

The experimental determination of the dynamic scat- 
tering law in liquid c esium is dated to th e early nineties 
l)Bodensteiner et all Il992t iGlaserl Il99l|) . Due to the 
relatively small incoherent cross section and to the low 
sound velocity value, after Rubidium liquid Cs is the 
most favorable alkali metal aiming at the study of col- 
lective dynamics by means of INS. Despite its late out- 
let, the work of Glaser and Bodensteiner reports an im- 
pressive state-of-the-art triple axis experiment, and a 



robust data reduction perf ormed with innovative algo- 
rithms l)Bodensteinerl . fl990 ) . 
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FIG. 24 Dynamic structure factor of Cesium at the melting 
point (circles). The continuous line i s the viscoelastic approx- 
imation. From ( Bo densteiner et ali 1992) 

The experiment is focused on the determination of 
the collective properties and, though a careful subtrac- 
tion of the incoherent contribution is performed, once 
more it emerges the intrinsic difficulty of determining 
the quasielastic part of the coherent spectrum. This 
notwithstanding, the data clearly show the departure of 
the collective dynamics from the strict hydrodynamic re- 
gion and the evolution toward the free streaming limit. 
This effect is quantified by the behavior of the FWHM 
of the quasielastic line reported in Figl25l while the hy- 
drodynamic prediction, based on purely adiabatic ther- 
mal fluctuations, predicts oJ\/2 = DtQ 2 , the actual 
linewidth is always below this limit in the whole ex- 
plored region, indicating the dominant presence of vis- 
cous processes in the quasielastic spectrum, as recently 
pointed out in other alkali m etals and liquid aluminium 
ijScopigno and Ruoccol l2004t) against an opposite inter- 
pretation in terms of l ineariz ed hy drodynamic models 
l|Singh and Tankeshwarl 120031 I2004D. The De Gcnncs 
narrowing was also observed ijDe Genneslll959|) . and the 
Q dependence of the linewidth was described within the 
hard sphere extended mode approximation fEa !671 sec- 
tion ^^3. The free streaming limit is not yet attained 
at wavevectors as large as twice the position of the main 
peak of the static structure factor ((Bodensteiner et~al\ . 
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1992). The lineshape analysis was performed with several 
approaches, within the e xtend ed hydrodynamic model 
llde Schepper and CohenL Il98(lf) . with the viscoelastic 
model (|Lovesevt Il97ll) and with two relaxation times ac- 
counting for both viscous and thermal processes. In this 
latter case, it was found a negligible role of the ther- 
mal process on approaching the first maximum of the 
structure factor, though the fitted values of the ther- 
mal relaxation time were in significant disagreem ent with 
the e xpected values 1/D T Q 2 (fig 9 of ref. l)Glaseit 
1991)). It was then pointed out the impossibility of 
discriminating the different models, given the s/n ra- 
tio of the available data. The Q— dependence of the 
longitudinal viscosity, extracted from S(Q,u> = 0) value 
(according to the prescription of generalized hydrody- 
namics), showed a decreasing be havior previously ob- 
served in Lennard - Jones systems ijAilawadi et all Il97ll 
Tankcshwar, 1994) and more recently in several other liq - 
uid metals ijMonaco et all l2004t IScopigno et al] . l2002b|) . 
Finally, the maxima of current correlation spectra showed 
the usual positive dispersion effect, which was interpreted 
as precursor of the solidification according to the usual 
idea of a transition between a liquid and solid-like regime. 
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FIG. 25 Linewidth of the coherent quasielastic spectral com- 
ponent in molten Cesium (circles) . The continuous line is the 
thermal value DtQ 2 , while the dashed line is the hard sph ere 
prediction from Eg. 1671 . From (Bodcnsteiner et al, 1992) 

Neutron scattering data on liquid cesium have also 
been used as a benchmark to develop an approach based 
on the id ea of timescale invaria nce of the relaxation 
processes llYulmetvev et aZ.ll200lf) . Within the Zwanzig 
Mori projectors formalism, one can construct an infinite, 
non Markovian, set of interconnected kinetic equations 
relatin g each mem ory function with the one of higher 
order l)Morilll965|) : 



dF(Q,t) 
dt 

dMi(Q,t) 
dt 

dMj(Q,t) 

dt 



-flj f drM\ (Q, t)F(Q, t — t) 
Jo 

Jo 

-Cl 2 i+1 f dTM i+1 (Q,T)Mi{Q,t 
Jo 



where F(Q, t) is the normalized density correlation func- 
tion and fli are characteristic frequencies of the process. 
Following the Bogoliugov approach of the reduced de- 
scription, one hypotizes the time scale invariance of the 
relaxation processes beyond a certain level, defining a clo- 
sure level Mj+i(£) w Mj(t) and thus getting an explicit 
expression for the DSF in terms of the spectral moments. 

To our knowledge no IXS experiments on liquid Cs 
have been reported, for the same kind of difficulties as in 
liquid Rb. 



B. Alkaline earth elements 

1. Magnesium 

Magnesium is one of the simplest divalent ele- 
ments. The coherent dynamic structure factor has 
been recently de t ermin ed at the SPring8 IXS beamlinc 
l)Kawakita et all [2003). The dispersion curve shows a 
8% deviation from the adiabatic sound velocity, with a 
maximum value lying halfway the hydrodynamic and the 
Coo value. An average relaxation time was determined 
(t = 0.094 ps), which is about one thir d of the one of 
the n eighboring alkali element liquid Na ijPilgrim et all 
1999). The analysis of the quasielastic line revealed a 
Q 2 broadening in the quasi-hydrodynamic regime, while 
around the De Gennes narrowing region the linewidth 
was success fully reproduced by the de Schepper and Co- 
hen model (|de Schepper et al\ . Il984[) . i.e. through Eq. 
(|67[) . Molecular dynamics simulations have been re- 
cently performed in this system by both classical molec- 
ular dynamics and orbital free ah initio simulations 
ijAlemanv et all . Il997t iGonzalez et al. i 12001ft . The two 
approaches give very similar results as far as the phonon 
dispersion is concerned, while the quasielastic contribu- 
tion is less pronounced in the ab initio calculation. In 
this respect, a compar ison with the experimental data 
ijKawakita et all 12003ft would be extremely interesting, 
though one should first convolute the calculated S(Q,u>) 
with the instrumental IXS resolution and take into ac- 
count for the quantum correction arising from detailed 
balance condition. 



C. Group III elements 

1. Aluminium 

One of the most puzzling results of early neutron spec- 
troscopy is the striking similarity between the spectra of 
polycrystalline and liquid Aluminium observed in Stock- 
holm in 1959 and publishe d in the final form a few years 
later IjLarsson et alJ . ll965|) . In this experiment TOF data 
were collected on a cold neutron spectrometer, but at 
those times multiple scattering corrections were almost 
impossible and therefore the spectra did not obtain a de- 
tailed explanation. A second TOF experiment was per- 
formed in the same period, but again the results were lak- 
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ing a detailed interpretation ((Cockingt Il967bj) . Fifteen 
years later the original experiment of Larsson was revis- 
ited, in an effort to re-analyze the results in the light of up 
to date theoretical developments. More specifically, the 
experimental re sults were used to test the old convolution 
approximation lIVinevardL Il9581. mean fi eld approaches 
llSingwi et aR 119701) and kin etic theory (jSiogreii Il978t 
ISioeren and Siolanderl Il978|) . Nevertheless, these data 
have been collected at very large energy and wavevectors 
value and always presented as time of flight scan, and 
therefore they are not very helpful to establish the ex- 
istence of collective modes in this system. Aluminium, 
indeed, is an almost purely coherent scatterer, but the 
high sound velocity value prevents the study of acoustic 
modes by means of INS (see table HJ. This is testified by 
a more recent IN S experiment performed at IN4 (ILL) 
ijEder et q/.l . ll98fj(l . Although multiple scattering correc- 
tion and constant -Q cuts reduction have been performed 
in this case, no evidence of collective modes could be re- 
ported due to the restricted kinematic region correspond- 
ing to the slow neutrons utilized (sa 55 meV). 
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FIG. 26 Dynamic structure factor of molten aluminium. 
Comp arison between the O F-AIMD calculation (continuous 
line) (|Gonzalez et all I20021 and the experimental IXS val- 
ues with their fitting based o n generalized hydrody namics 
iScotwno et »i.t 12001 al). From JOonza.lez et aZlEimf) 



Much of the knowledge about the microscopic dynam- 



ics in liqu id Al relies , in fa ct, on the numerical work of 
Ebbsjo ijEbbsio et a?l ll980h who calculated the dynamic 
structure factor utilizing two different local pseudopoten- 
tials and the local pseudopot ential originally developed 
by Ashcroft l)Ashcroftl Il966|) . The dynamic structure 
factor has been shown to be somehow reminiscent of the 
viscoelastic model with the addition of a gaussian term, 
accounting for the approach to the high Q, free streaming 
limit. He was then able to predict the existence of collec- 
tive modes for Q < 10 mif 1 , though he reported sound 
velocity values distinctly larger then the adiabatic value 
over the whole explored range (Q > 3 nm _1 ). Triggered 
by this observation, a modifi ed version of th e viscoelas- 
tic approach was developed <|C T a,skelH |l986) and tested 
against the data of Ebbsjo. At the same time, it has been 
proposed a single relaxation process scenario based on a 
sech m emory function shape wi thin the Mori-Zwanzig 
scheme l|Tankeshwar et a/lll988j) . 




FIG. 27 Sound dispersion of liquid aluminium from the max- 
ima of the curren t correlation function : open circles, OF- 
AIMD calculation llGonzalez et al. , 12001). full circles IXS ex- 
perimental values IScrjpigno^^Z , 12001a) . The dispersion 
from the maxima of the dynamic structure factor numerically 
evaluated is also reported (open triangles), a s well as the hy- 
drody namic value (continuous line). From iGonzalez et all 
2002) 

The first experimental observation of collective modes 
in liquid Al ha s been reported much m ore recently by 
means of IXS llScopigno et all l2001a|) . An high fre- 
quency regime has been observed for Q > 5 nm _1 , while 
below this value the dynamics approaches the hydrody- 
namic limit, though the transition is not fully accom- 
plished at the lowest investigated wavevector, Q = 1 
nm _1 . The scenario arising from the IXS study is much 
similar to the one characterizing alkali metals, though 
with significant quantitative differences such as a more 
intense quasielasic component testifying a more impor- 
tant role of the structural relaxation in this system. 

The IXS data have been recently use d to test orbital 
free ab initio calculations (OF-AIMD) ((Gonzalez et all 
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FIG. 28 Resolution deconvoluted, classical li neshape ut ilized 
to described t he IXS spectra of molten Li llScopigno et all 
l2000d . l2001al . l2002aft . Na and Al, reported on relative mo- 
mentum and energy scale (see text) 

1200 lL l2002t) . The overall agreement is quite satisfac- 
tory though the numerical calculations show somehow 
lower sound velocity value and tend to overemphasize 
the inelastic components. From FigE|one can argue the 
importance of the dynamical variable representing the 
sound velocity. The presence of positive dispersion, in 
particular, is strongly affected by the choice of the max- 
ima of Cl(Q, rather than the ones of S(Q, u>). 

In Fig. [2H1 finally, we report a comparison of the line- 
shape obtained from the resolution-deconvoluted, clas- 
sical lineshape utilized to fit the IXS spectra of liquid 
lithium, sodium and aluminium at the same reduced val- 
ues of momentum (Q/Qm) and energy (toto — uj/loq 
with to = \/m/kBT m /QM) transfer. As can be clearly 
evinced, the attitude of alkali metals to sustain density 
fluctuation is much more pronounced than in other sim- 
ple liquid metals. 



2. Gallium 

Among simple liquid metals, Ga is endowed with pe- 
culiar structural and electronic properties. In addition 
to the low melting point (T m = 303 K), it shows an ex- 
tended polymorphism in the solid phase with complex 
crystal structures where a competition between metallic 
and covalent bindings takes place. Despite the nearly free 
electron electronic DOS anomalies associated with some 
covalency residue are present. Moreover, the first peak 
of the S(q) presents a hu mp characteristics of non clo se- 
packed liquid structures l)Bellissent-Funel et all \l989). 

Early inelastic scattering on liquid Ga were per- 
formed at the beginning of the seventies with neutrons 
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llBosio et oil Il97fit lOlaser et all Il973t iLoffleit Il973l 
iPage et allll973fl . Due to the quite large sound velocity 
value, compared to the available kinematic range, these 
studies were mainly addressed to the investigation of the 
quasielastic part of the dynamic structure factor. 

Twenty years later, another series of INS exper- 
iments were perform ed on a triple axis instrument 
l)Bermeio et all Il994f) with the aim of ascertaining the 
presence of low Q collective modes just above the melt- 
ing temperature. Indeed, although by virtue of the above 
mentioned anomalies liquid Ga seems to elude the pic- 
ture of the high frequency dynamics emerging in all the 
monatomic liquids, on the basis of the shape of the inter- 
action potential evidence for collective modes should be 
expected below Q m /2. Quite surprisingly, no evidence of 
inelastic signal was reported in the region were longitu- 
dinal modes were expected on the basis of the hydrody- 
namic sound velocity. This result was interpreted as an 
overdamping effect traced back to the high value of lon- 
gitudinal viscosity (« 10 cP). Additionally, by compar- 
ing constant energy scan (see Fig l29|) to the expression 
of Buchcnau for acoustic-like plan e wave excitations in 
amorphous solids ijBuchenaul 1 1 985j) . an excess of scatter- 
ing was reported for frequency distinctly larger then the 
maximum frequency of the acoustic branch. This result 
was interpreted as the fingerprint of the presence of high 
energy optic-like modes. 

A few years later, then, a new set of experiments were 
performed at h i gher t emperatures by the same authors 
[|Bermeio eTal\ . Il997|) . and a contrasting behavior with 
the previous findings was reported. More specifically, 
the appearance of non-overdamped sound modes was re- 
ported, accompanied by a second, higher frequency mode 
of presumably optical origin. The discrepancy between 
the low and high temperature experiment was ascribed to 
a viscosity drop of a factor « 7 and therefore to a narrow- 
ing of the acoustic mode according to the hydrodynamic 
expression II 71 

The presence of an higher frequency mode appearing 
at in the constant Q scan for wavevectors larger then Q m 
seems to corroborate the presence of an optic-like exci- 
tation sugges ted by the previously mentioned constant 
energy scans ijBermeio et all[l99^ . 

A recent IXS experiment performed on liquid Ga llium 
just above the melting point (Sc opigno et oiJ . l2002bl) por- 
traits the collective dynamics in a much similar fashion 
to the one of Alkali metals and of liquid Al, thus re- 
moving the anomaly suggested by the neutron experi- 
ments. Collective modes, in fact, have been unambigu- 
ously observed in the low temperature region where neu- 
trons suggested an overdamped regime. This result sug- 
gests the inadequacy of the expression 1171 to estimate 
Brillouin linewidths, which can be easily understood in 
terms of the gener alized hydrodynamics re sults reported 
for alkali metals l|Scopigno et all l2000b|) : outside the 
truly hydrodynamic region, the viscous relaxation dy- 
namics proceeds over two distinct physical mechanisms, 
the structural relaxation and the short-lived rattling dy- 
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FIG. 29 INS constant energy scans (circles) compared to 
the model of Bu chenau for plane wave excita tions in solids 
l|Buchenaii 1198.1 ) . From ijBermeio et al.Lll994l) . 



namics. On the high frequency region of interest, the 
structural relaxation is frozen (the system is respond- 
ing as a solid) and therefore the viscosity associated to 
this process does not contribute to the sound damping. 
Moreover, the thermal contribution in Eq. I|17|) might not 
be correct at wavevectors as large as a few nm _1 , since 
the adiabatic regime could be r eplaced by as isother- 
mal o ne, as already pointed out ijScopigno and Ruoccol 
l2004j) . Consequently, Eq^lis an overestimate of the ac- 
tual linewidth (which in the c ase of liquid lithium h as 
been quantified as a factor two ijScopigno et all 12 000b)). 
In Fig. 1301 we report the two experiments for simi- 
lar values of t he momentum transf e r. The viscoe las- 
tic prediction l)Bermeio et all 11994 lLovesevt Il97l|) is 
also reported, showing how it clearly underestimates the 
quasielastic contribution, though it provides a reasonable 
estimate of the genuine inelastic mode. The IXS findings 
have been recently corroborated by a new accurate INS 
experiment, used to test the reliability of a model in- 
ter action_p^tarti^bycomparing the dynamic structure 
factors l)Bove et all l2005|) . 

Summing up, the lack of low temperature collective ex- 
citations reported in this system with neutrons is prob- 
ably due to the difficulty of a reliable extraction of the 



FIG. 30 Compariso n between the INS (stars ijBermeio et all 
1 19941) and triangl es (iBove et all 20051)'! a nd IXS (open circles 
with error bars <|Scopigno et all 1200 2b)) determinations of 
S(Q,oj) in gallium at the melting temperature, for two differ- 
ent values of the (constant) mom entum transfer. The dotted 
line is the viscoelastic prediction iBermeio et aZlll994f) (con- 
voluted with the INS instrumental resolution and accounting 
for the detailed balance condition), while the continuous line 
is the best fit to the IXS data utilizing a memory function ac- 
counting for the thermal relaxation and two viscous processes 
(see text). Molecular dynamics symulations for the coherent 
(dotted line) and total (dash-dotted line) S(Q,cu), convolut ed 
to the INS resolution, are also reported iBove et all 120051 . 



coherent part of the scattering function. On the other 
side, the interesting observation of possibly optic-like 
high frequency modes certainly calls for further inves- 
tigation although, in our opinion, does not justify per se 
a description of the acoustic dynamics in terms of crys- 
talline reminiscent dynamics. The analogy with several 
IXS results on different systems suggests indeed a promi- 
nent role of the topological disorder in characterizing the 
acoustic branch. 

Very recently, an IXS experiment explored the high Q 
region (20 < Q < 70 nm" 1 , i.e. lengthsc ales smaller than 
the s ize of the first coordination shell ijScopigno et all 
120051) . While generalized hydrodynamics provides a co- 
herent picture of the dynamics in the lower Q region, not 
much is known about collective dynamics at such short 
lenghtscales. For hard spheres, Enskog's kinetic theory 
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FIG. 31 Full width at half maximum of the dynamic struc- 
ture factor as determined by IXS in the kinetic regime (open 
circles). The prediction according to Enskog's theory is also 
shown (continuous line) (IScopigno et all 120051) . 
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predicts in this region the dominant effect of a general- 
ized heat mode. Liquid Gallium, however, by no means 
can be modelled as an hard sphere fluid, for the above 
mentioned structural and electronic properties. Surpris- 
ingly, it turned out that a description in terms of heat 
mode (Eq. still applies, at the price of introduc- 

ing an effective hard sphere diameter (larger than the 
one associated to the first S(Q) maximum), which prob- 
ably accounts for the associative tendency of this liquid 
(dimer-like structures). 

D. Group IV elements 

1. Silicon 

Due to its several unusual properties, liquid Si is al- 
ways classified as a non simple liquid. While in the 
crystalline phase Si is a diamond structure semiconduc- 
tor, it undergoes a semiconductor-metal transition upon 
melting, which is accompanied by a density increases of 
about 10%, and by significant structural changes (the 
coordination number grows from four in the solid state 
to about seven in the liquid). Similarly to Gallium, the 
static structure factor, S(Q), exhibits a shoulder on th e 
high-Q side of the first peak ijWaseda and Suzuki Il975|) . 
a feature that cannot be reproduced using a simple hard- 
sphere model, appropriate for alkali metals. 

No neutron scattering data exist to the best of 
our knowledge, while very recently two Inelastic X 
ray Scattering experiments have been p erformed both 
at the ESRF jHosokawa et all l2003a|) and SPring8 
l|Hosokawa et a/.tl2003d) . 

In the first of the two above referenced experiments, 
a positive dispersion of 15% has been found. In the 
second experiment, an higher resolution and more ac- 
curate study has been carried out, allowing to follow 
the transition from the high frequency to the low fre- 



FIG. 32 Sound velocit y and attenuation in molten silicon 
jHosokawa et all l2003cfl 



quency regimes. The data were analyzed in terms of 
DHO model for the inelastic component and Lorenzian 
for the quasielastic, and no significant quantitative differ- 
ences were detected utilizing the same memory function 
scheme proposed for ot her IXS studies on liquid metals 
l)Sconigno et all l2000b|) . 

In the vicinity of the main peak of the static structure 
factor, the lorenzian shape for the quasielastic component 
turned out to be inadequate, and has been replaced by 
a combination of Lorenzian and gaussian contributions 
(pseudo-Voigt). 



2. Germanium 

Liquid Ge shares the same peculiarities as liquid Sili- 
con, though with some slight quantitative differences. 

IXS data on liquid Ge have been recently obtained 
l(Hosokawa et a/1 l200ll) . and they show evidence for col- 
lective propagating modes. An analysis based on a Loren- 
zian shape for the quasi-elastic and a DHO for the in- 
elastic modes, revealed the absence of positive dispersion 
effects in the investigated Q range (2 28 nm" 1 ). On 
our opinion, this result calls for further investigations, as 
this is an almost unique feature in respect to the other 
monoatomic liquid metals investigated sofar (especially 
Silicon, which has very similar structural properties). 
The De Gennes narrowing has been analyzed in terms 
of extended hydrodynamic heat mode, utilizing the an- 
alytical e xpression obtained within hard-sphere approx- 
imation l)Cohen et all Il987j) . but the quite large error 
bars and the limited spanned Q range prevents to draw 
a final conclusion. 



40 



3. Tin 

Tin is the heaviest of 4B element. Its structural prop- 
erties are quite similar to the one of alkali metals, with 
a coordination number close to 12 but with the typi- 
cal shoulder on the high Q side of the main S(Q) peak 
ijWaseda and Suzuki ll9T5^ which is typical of Si, Ge and 
Ga. 

The first INS experiments in liquid Tin have to be 
traced back to the early sixties. Similarly to other early 
neutron scattering experiment no clear picture of the 
microscop i c dyn amics could be outlined. In some case 
llPalevskvl Il96lft the vineyard approximation ijVinevardl 
Il958ft was used to analyze the data, while in another 
study experimental str ategy for suppres s ing m ultiple 
scattering were tested l|Brockhouse et all ^963). Al- 
though a w avevector-energy plot was obta ined from raw 
TOF data l|Cocking and Egelstaft ll965tJ) . according to 
Copley and Lovesey no side peaks should be observed 
trasforming TOF data to S(Q,ui) on constant Q slices 
ijCoolev and Lovesevlll975h . 

Constant Q IXS spectra of liquid tin have been very 
recently reported for low (T = 593K) and high (T = 
1273-fQ temperatures at the ESRF l)Hosokawa et all 
2003b!). The sound velocity seems to exceed the hydro- 
dynamic value at both temperatures of 6% andl2%, re- 
spectively. This notwithstanding, these quantitative es- 
timates must be taken with care, due to the quite large 
error bars. Moreover, the dispersion curves have been de- 
termined from the DHO frequency, neglecting the effect 
of the quasielastic component. For Q values close to the 
main peak of the struc t ure fac tor, as in the case of liq- 
uid Si ijHosokawa et a/.l . l2003cft . the lineshape of S(Q,ui) 
turned out to be empirically described by a combina- 
tion of Lorentian and Gaussian contributions or, equiva- 
lently, by a memory function analysis similar to the one 
reported for alkali metals (?). As a matter of fact, at such 
large Q values, as previo usly observed for liquid lithium 
l)Scopigno et aZJ . l2000albl) . the microscopic dynamics un- 
dergoes a transition from the (generalized) hydrodynamic 
behavior to the free streaming limit and a detailed de- 
scription of such transition is still missing. 



4. Lead 

Molten lead has been one of the first metals to 
be investigated by Inelastic Neutron Scattering, as 
the first experiment can be traced back to the fifties 
l)Brockhousel Il955l lEeelstaflt Il954l). Details of the ex- 
periments performe d up to 1975 llCocking and Egelstafl , 
1968J iDorner et all Il965l Il967t iRandolph and Singwi 
1966^1 have been exhaustively rev iewed by Copley and 
Rowe l)Coplev and Lovesevt Il975ft . One of the most in- 
teresting results has been the evidence of both a longitu- 
dinal and a transverse branch in the dynamic structure 
factor, though this result was presented with some cau- 
tion as the transverse mode could be an artifact arising 



from multiple scattering effects. 

In the early eighties new INS studies performed 
with both TAS and TOF spectr o meter s appeared 
ijSoderstroml Il98ll ISoderstrom et all Il980ft . aiming to 
validate the presence of a dispersion relation and of a 
transverse branch. A longitudinal mode, compatible with 
the hig her frequency excita tion previously reported by 
Dorner l)Dorner et all Il967ft was reported, whose sound 
velocity is consistent with hydrodynamic value. The ac- 
cessible kinematic range is too limited to asses any ev- 
idence of positive dispersion effect. No evidence for a 
lower frequency excitation was instead reported, corrob- 
orating the hypothesis that such feature is an artifact 
stemming from multiple scattering. 

Liquid lead has been recently in focus of molec- 
ular dynamic simulations aiming to describe collec- 
tive dynamics in t e rms of generalized kinetic modes 
ijBrvk and Mrvgloaj. l2001albft . Beyond the hydrody- 
namic region, different branches corresponding to sound 
and heat waves have been identified, and their nature has 
been extensively discussed. 



E. Group V elements 

1. Bismuth 

The first inelastic scattering investigation of molten 
bismu th was originally reported by Cocking l)Cockind . 
Il967al) . who reported a dispersion curve extracted by 
TOF neutron spectra. At the IAEA symposium of 
1968 two experiments Bi were presented: in one case 
two dispersion branc hes were obtained from TOF data 
ijTunkelo et a/.l . 11968ft . though the low frequency excita- 
tion was probably due to an artefact of a missing multiple 
scattering correction. In the s econd study data were con- 
verted to constant Q spectra l)Mateescu et a/.Lli968ft . 

Liquid Bi h as been recently the subject of new INS 
investigations l)Dahlborg and Olssonl ll 982ft . Measure- 
ments were performed just above melting at T = 578 
K, in a wavevector region spanning from just below the 
first maximium in S(Q) up to wavevector as high as 
70 nm _1 . This kinematic region lies above the region 
were collective modes could be expected, so the study 
mainly deals with the quasielastic spectral component. 
Generalized hydrodynamics models based upon a sin- 
gle relaxation time were tested against the experimen- 
tal data utilizing differ ent approximation for the mem - 
ory function shape ij Ailawadi et a/.L ll97l1: lLovese 3 Ell- 
Experimental data were then unfolded by instrumental 
resolution modelling the quasielastic shape as Lorenzian 
and the resolution as gaussian. The main point of the 
paper is the determination of the Q dependence of the 
spectral FWHM, which is also compared to other sys- 
tems. In particular, it is pointed out how, beside the 
expected De Gennes narrowing occurring at Q^ni the 
FWHM shows a minimum rather then the expected max- 
imum at Q ~ l.bQm (which characterize the dynamics 
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FIG. 33 FWHM of molten Bi compared to other simple liq- 
uids. It can be noted the presence of an intermediate min- 
imum between the first two maxima which is not present in 
the other simple fluids. From iDahlborg and Olssonl Il982f) . 

in several simple liquids such as Pb, Rb and Ar). This 
anomaly is related by the authors to the shoulder ob- 
served in the static structure factor just above Q m , and 
therefore identified as a non-simple nature of liquid Bi. 

The kinetic region has been investigated in liquid 
Bi within the generalize d collective mode approach 
l)Brvk and Mrvglodl l2000|) . The presence of high fre- 
quency kinetic branches has been ascertained, and it has 
been pointed out that their weight is too small to make 
them visible in the dynamic structure factor. This result 
seems to be in agreement with recent IXS findings on 
a very similar system, namely liquid Ga, in which on ly 
acoustic modes were reported llScopigno et aZJ . l2002b^ . 

F. Transition metals 



of 5E sa 1 meV. The data are analyzed with an empirical 
model consisting of a DHO for the purely inelastic part 
and the sum of two lorentian functions accounting for the 
quasielastic contribution. While the inelastic component 
is no doubt of coherent nature, the narrower of the two 
lorentians is ascribed to incoherent scattering, and mod- 
elled as a simple diffusive term of linewidth DQ 2 . Given 
the values of D, this results in a quasi elastic contribution 
which could not be resolved by the much broader reso- 
lution function. The linewidth of the broader lorentian 
is almost Q independent, and its origin is ascribed by 
the authors to an incoherent process, on the basis of the 
coincidence of the DHO area with independent (static) 
structure factor determinations. Turning our attention to 
the collective dynamics, the extrapolated high frequency 
value of the sound velocity (coo(Q — ► 0) = 2100±80 m/s) 
, obtained by the DHO frequency parameter (therefore 
approximately equal to the maximum of the current cor- 
relation function, the difference being due to potential 
quasielastic coherent contribution), turns out to largely 
exceed the hydrodynamic value (c s = 1451 m/s), suggest- 
ing a huge positive dispersion effect close to 50% largely 
exceeding similar effects reported in other metals. This 
result is rationalized in terms of Bohm-Staver model, 
which provides the expression of Eq.J7SJ) for the sound 
velocity yielding, for molten Hg, ci(Q — * 0) » 2090. 
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1. Mercury 

Experimental studies of microscopic dynamics in liq- 
uid mercury are very recent compared to the systems 
reviewed so far, and they were presented at the LAM XI 
conference (Yokohama, Japan). The first investigation 
was obtained by means o f INS at the INI facility of the 
ILL dBove et all . l2002albD . In this work, an detailed in- 
vestigation of the dynamic structure factor is undertaken 
at room temperature, and presented as constant Q cuts 
in the range 2.5-^12 nm _1 with a high energy resolution 



FIG. 34 DSF of molten Hg measured by means of INS at the 
indicated Q values. The inelastic and the quasielastic compo- 
nens, modeled with t wo lorenzians and a DHO, respectively, 
are also shown. From jBove et all l2002al) 

Nearly in coincidence with the neutron experiment of 
Bove et al., an IXS study of the coheren t dynamics in 
liquid Hg appeared ijHosokawa et a?J . l2002|) . in which the 
dynamic structure factor is examined in a wavevector re- 
gion extended up to Q rn ~ 25 nm , with a factor 2 
worse resolution, but in a significantly larger energy re- 
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gion. As in the work of Bove et ai, the genuine inelastic 
features of the data are modelled with a DHO function, 
but it appears clearly from the raw data that a coherent 
quasielastic term dominates the u> rs region. This lat- 
ter contribution is modelled with a lorenzian shape. Al- 
though the authors neglect the presence of this quasielas- 
tic in the calculation of the sound velocity value, taking 
the DHO frequency as the relevant parameter they ob- 
tain a value of c^Q — > 0) = 1840 m/s, an estimate 
which is directly comparable with the corresponding INS 
determination. This discrepancy, which may be due to 
the limited energy range at the low Q's of the INS exper- 
iment, as well as to the non negligible resolution effect on 
the lowest Q's of the IXS experiment, calls in our opin- 
ion for further investigations and suggests to take with 
some care any interpretations in terms of Bohm Staver 
model of the positive dispersion at least by the quantita- 
tive point ov view. Hosokawa et ai, on the other hand, 
cast the anomalous dispersion in the usual framework 
of the shear relaxation. The IXS work also confirms the 
non negligible presence of quasielastic signal in the coher- 
ent dynamic structure factor, suggesting that it should 
be taken into account also in INS data treatment. We 
expect that the broader lorenzian contribution reported 
by Bove et al., for example, could be at least partially 
coherent in nature. 
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FIG. 35 Left panel: IXS determination of the DSF in Hg 
near the melting temperature. Right upper panel: dispersion 
relation and sound attenuation properties as deduced by a 
DHO + one Lorenzian fit. The low and high frequency limits 
are also reported . Right lower panel: Corresponding sound 
velocities. From llHosokawa et aZ.Ll2002l) . 

A more recent INS study contributed to shed some 
light on the possible origin of th e quasielastic spectral 
components IjBadval et all |2003|) . The experiment has 
been performed on the TOF spectrometer MARI, op- 
timized to access a restricted kinematic region (—6 < 
lj < 6 meV at the lowest accessed momentum trans- 
fer Q ss 4 nm" 1 ) with an energy resolution sufficient to 
study the diffusive dynamics (SE = 0.4 and 5 = 0.8 meV 



at the two incident energy utilized). The narrower in- 
coherent contribution (self diffusion), resolution limited 
in the experiment of Bove et al., was now detected with 
a pro cedure similar t o the one applied in liquid potas- 
sium (fCabrillo et all 120021) . i.e. fitting the data with 
two lorenzian in a restricted wavevector region where the 
two quasielastic features are well separated, determin- 
ing the analytical Q dependence of the FWHM of the 
diffusive term, and finally focusing on the broader com- 
ponent over the whole momentum transfer region keep- 
ing all the parameters of the diffusive term fixed. With 
this procedure, the narrower lorenzian is confirmed to be 
incoherent in nature, and well described by Fick's law 
prop erly modified accord i ng to the revised Enskog's the- 
ory ijKamgar-Parsi et all Il987|) . As far as the broader 
component is concerned, the authors point out that more 
than one lorenzian is needed to describe it at increasing 
wavevector, then they discuss its possible origin. First, 
they point out that the thermal origin of this broad com- 
ponent is ruled out by its low Q intensity, which is by 
far larger than the one expected by the Landau Plazeck 
ratio. Moreover, the linewidth reaches a constant value 
on decreasing wavevector, instead of following the Q 2 de- 
pendency of the heat mode. Second, they examine the 
possibility of a cag e diffusion mecha nism, as proposed 
in MD simulations ijBove et aZJ . 12002^1 . In this respect, 
they point out how the experimentally observed mode 
intensity is too large than expected, but they propose an 
enhancement mechanism based on valence fluctuations 
which could be active at low wavevectors amplifying the 
expected intensity. 

Very recently, new state-of-the-art IXS experiments 
have been reported in expanded mercury near the criti- 
cal point (T c — 1751 K, P c = 1673 bars and p c — 5.8 g 
cm -3 ), aiming at the investigation of collective dynamic s 
at the metal-non metal transition ijlshikawa et all 12004). 
Despite extremely difficult experimental conditions, the 
speed of sound has been accurately measured and no 
significant changes have been observed in the transition 
from the metallic (mp c = 13.6 g cm~ 3 ) to the insulating 
(mp c = 9.0 g cm~ 3 ) phase, while the ultrasonic sound ve- 
locity exhibits a significant drop across the same thermo- 
dynamic point fYao et aiJ .IT996). Only upon further ex- 
pansion in the insulating phase the high frequency sound 
velocity ultimately drops reaching the adiabatic value. 
This result seems to indicate the presence of very large 
positive dispersion as peculiar of the metal-non metal 
transition, opening a new experimental route to the in- 
vestigation of the interplay between acoustic and trans- 
port properties. 

Summing up, though the microscopic dynamics in 
molten Hg has been the subject of deep investigations 
in the last few years, some aspects still remain contro- 
versial. The sound velocity as determined by INS and 
IXS are not fully consistent with each other, though both 
techniques show positive dispersion effect the INS result 
show a very large effect never observed sofar. On the one 
side it has been emphasized the role of electronic proper- 
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ties (jBove et a/.ll2002a|) . while, on the other side, the col- 
lective dynamics as determined by IXS closely resemble 
the o ne of several other simple fluids ijHosokawa et all 
2002). The most intriguing aspect concerns, how- 
ever, the interpretation of the quasielastic component 
of S(Q,ui). Neutron scattering d ata suggest a negligi - 
ble effect of thermal fluctuations l|Badval et all [2003), 
addin g a piece of information to a recently debated is- 
sue llScopigno and Ruoccol 12004 ISingh and Tankeshwarl 
l2QQal2004). On the other side, the IXS data unambigu- 
ously show the presence of a coherent quasielastic dynam- 
ics, which no doubt has to show up in neutron experiment 
as well. The broad quasielastic component as observed 
with two different experiments s how however opposi te Q 
trends, mo notonically inc r easing ijBove et a/.Ll2002"al) and 
decreasing ijBadval et a/.l 12003]) . respectively. Badyal et 
al. suggest a cage diffusion process, enabled by valence 
fluctuations mechanism. We believe that the broad mode 
observed in INS could be closely related to the coherent 
quasielastic scattering reported in the IXS data. In this 
case, the similarity with several other investigated sys- 
tems would suggest an interpretation in terms of a high 
frequency structural relaxat i on process llMonaco et a/1 
120041: IScopjgno et ~aH l2000d l2001al l2002alblh Very re- 
cent IXS investigations, however, suggest an enhance- 
ment of positive dispersion at the metal non metal tran- 
sition, pointing out how changes in electronic trans- 
port properties dram atically affect acoustic properties 
illslnkmvn ,1 aILl200l. 



2. Nickel 

Early investigations in liqui d Nickel have been re ported 
in 1977 with TOF technique l|,Iohnson et n/.lll977l) . Two 
different isotopic concentrations, one with the natural 
abundance ratio and the other a wholly incoherent mix- 
ture, were chosen in order to study separately the coher- 
ent and incoherent scattering contributions. These lat- 
ter, in turns, have been related each oth er through the 
Vineyard approximation ijVinevardl Il958) . The spectral 
density g(uj) was then extracted from the low Q limit of 
the self dynamic law. The coherent S(Q,oj = 0) were 
reported for wavevectors above Q — 22 nm" 1 , therefore 
beyond the first brillouin pseudozone. Consequently, this 
study was not able to ascertain the existence of low Q 
collective modes. 

A neutron scattering experiment has be en recently per- 
forme d at the INI facility of the ILL (|Bermeio et a,l\ . 
.2000). Constant Q spectra have been collected from 
Q = 8 on" 1 all the way up to wavevectors as high as 
twice the main peak in the static structure factor. The 
data have been analyzed utilizing two lorentian terms for 
the quasielastic (coherent and incoherent) contribution 
and one DHO function for the purely inelastic spectral 
features. Evidence of collective modes has been reported 
in the whole investigated Q domain, despite the relatively 
high viscosity value (rj s 5.7 mPa s, i.e. one order of 



magnitude larger than the typical values for alkali met- 
als). The reason for this apparent oddness can be traced 
back to arguments similar to the one applying to the case 
of gallium, i.e. to the freezing of the diffusional motion 
over the probed high frequency regime, which reduces 
the effective Brillouin damping in respect to the hydro- 
dynamic prediction of Eq. I|17ll . 




FIG. 36 Sound velocity in liquid Nickel determined by INS 
(full circles) and MD (open circles). The isothermal (ct) and 
adiabatic (c s ) values are also indicated. 

More interestingly, the low Q limit of the sound veloc- 
ity seems to approach the isothermal and not the adia- 
batic value, as shown in Fig. 1361 It is worth to stress 
how, given the large value of 7 (and therefore the large 
differences between ct and c s ), this observation strikingly 
holds beyond the reported experimental error. If con- 
firmed, this result would substantiate the hypothesis of 
an intermediate isothermal domain bridging the hydro- 
dynamic limit and the high frequ ency regime, which has 
been recently matter of debate llSconigno and Ruoccol 
120041: ISingh and Tankeshwarl 120031 l2004h . On the other 
side, previous molecular dynamics simulations indicate 
higher values of the sound v elocity which agree qu ite well 
with the adiabatic response <)Alemanv et a/.lll9 98). Sum- 
ming up, liquid Ni seems to be an ideal system to test 
the evolution of the thermal relaxation once the hydrody- 
namic limit is abandoned. An IXS investigation would be 
helpful to clarify this issue, though the small absorption 
length, the high melting temperature and the reactivity 
of Ni poses severe experimental challenges. 



3. Copper 

Time of flight neutron scattering data have been re- 
ported for this system long time ago in solid and liq- 
uid phase. The accessed kinematic range was quite 
broad (Q > \Qnm~ 1 and E < 30 meV) and the 
data have been analyzed w ithin pioneering models 
iEgelstaff and Schofieldl Il962|h More recently, an IXS 
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experiment has been performed at the ESRF, though this 
work is still in progress a very preliminary estimate of the 
sound velocity gives a value of 4230 ± 70m/s, well above 
the hydrodynamic value. 



G. Solutions of metals 

Alkali metals easily dissolve in water, molten alkali 
halides and ammonia, resulting in a free electron and a 
positively charged ion. Given the relatively low electronic 
concentration of the saturated solutions, these are ideal 
systems to challenge the validity extents of plasma-based 
theories introduced in section Hi .HI and the relative ap- 
proximations for the dielectric function. Models like the 
RPA, indeed, are expected to hold in systems with low 
r s (or, equivalently, high electronic density) such as bulk 
metals, while they reach their limits on increasing the r s 
value. 

A second interesting aspect concerns the presence in 
liquid metals of the so called Kohn anomaly, i.e. a kink in 
the dispersion curve occurring at Q — 2kp and reflecting 
a singularity of the dielectric function which is observed 
in metallic crystals, but which is not yet established 
in the disordered phase. The high frequency dynam- 
ics of metal-ammonia systems, in particular, have been 
recently investigated by means of X-rays and neutrons. 
High resolution IXS performed in lithium-ammonia solu- 
tions allowed to detect high frequenc y excitations, soften- 
ing at twice the Fermi wavevector kf ijBurns et alJ .120011. 
Unfortunately, at the investigated concentrations 2k p is 
close to Q — Qrm he. the main peak of the static struc- 
ture factor. This coincidence generates an ambiguity in 
the interpretation of the observed dip in the sound dis- 
persion, as it is not clear whether this feature is related 
to the structural periodicity as in non ordinary fluids or 
has to do with a purely electronic effect. The reported 
sound velocity values are intermediate between the one 
of the bare ions and the one appropriate for the Li(NH3)4 
vibrati ng network , thoug h close to this latter. In a later 
study l)Said et all l2003|) . collective excitations are ra- 
tionalized in terms of ionic plasma oscillations, and the 
sound velocity values measured at several metal concen- 
trations are compared with the prediction of the Bohm- 
Staver expression of Eq. 1751) , taking as relevant mass ei- 
ther the bare ionic value or the metal-ammonia unit. In 
both cases the predictions do not agree with the mea- 
sured values, and this seems to be a signature of the 
well known failure of the RPA approximation in the low- 
density regime, where the electron-electron interactions 
are releva nt. A contemporary I NS study on deuterated 
ammonia l|Sacchetti et all l2003|) also addresses the devi- 
ations from the BS model. In this case an improvement is 
achieved accounting for both the finite ionic size and the 
electron-electron interactions. As far as the first aspect 
is concerned, an alternative renormalization of the free 
ionic plasma frequency is undertaken, while an expres- 
sion for the dielectric screening function, going beyond 



the RPA approximation of Eq. (|74|l and accounting for lo- 
cal fields effects, is proposed. As far as the Kohn anomaly 
is concerned, this study suggests that the value of 2k f 
(w 10 nm -1 in saturated lithium- ammonia solutions) has 
to be compared with Q ss 20nm~ 1 , which is actually the 
second peak of S(Q), related to the N — N periodicity, 
rather than with the first peak (Q m « 20nm _1 ) related 
to the Li(NBi3)4 periodicity. This might arise from the 
different X ray and neutrons cross sections: while in the 
first case the two peaks have similar intensities, the neu- 
tron diffraction strongly enhance the N — N peak. A 
later study has shown how an unambiguous separation 
between 2Kp and Q m occurs in low density solutions, 
although in this case the kink observed i n the disper- 
sion curve is almost within the error bars l|Giura et all 
2003). Further studies, therefore, seem to be necessary 
to draw a conclusive picture about the presence of the 
Kohn anomaly in metallic fluids. 
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FIG. 37 Random phase approximati on at work: data taken 
from Table U and l|Bove et qlj |20 05|) for pure elements and 
from (iGiura et all 12003: lSeddefa/Tl200"3T) for alkali-ammonia 
solutions, with the BS values estimated through Eq. 1751 . 
Full circles: the measured sound velocity, c, is the maximum 
high frequency value determined over the Q range probed by 
inelastic scattering techniques, i.e. includes the positive dis- 
persion effect. Open circles: c is here the adiabatic value. The 
subscript of Li and Na indicates the concentrations in ammo- 
nia solutions. The Na group is relative to several different 
temperatures. 

The validity extent of the RPA approximation for the 
determination of the sound velocity via the Bohm-Staver 
expression (|75[1 is depicted in fig. 1371 in which we re- 
port the relative deviations of the BS calculated values 
(crpa) from the experimental ones, for systems of dif- 
ferent r s , ranging from pure metals to alkali-ammonia 
solutions. In the latter case (r s > 6) the BS predictions 
underestimate more than 50% the calculated values. The 
deviations, however, show a trend which monotonically 
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decreases towards low r s elements and finally get nega- 
tive for r s values close to 2. 

A final remark concerns the choice of the dynamical 
variable to calculate the sound velocity when one is look- 
ing at subtle effects as in the present case. First, ac- 
cording to the hydrodynamic definition of sound velocity 
in liquids, one should look at the maxima of the cur- 
rent correlation spectra. While the difference with the 
maxima of S{Q,w) is usually negligible in crystals, in- 
deed, there can be a significant discrepancy in strongly 
overdamped cases such as the one metallic solutions. In 
this respect, the choice of the DHO to describe S(Q,u>) 
implicitly overcame this problem, as the characteristic 
frequency of this model coincides indeed with the max- 
ima of Cl{Q,oj). Second, in all the reported studies, 
the "phonon" velocity is extracted through ad hoc mod- 
els (DHO, extended hydrodynamic model etc.) looking 
only at the genuine inelastic component. Again, in liq- 
uids, the full Cl(Q, uj) spectrum should be considered, as 
when S(Q,uj) has a significant quasielastic contribution 
this latter can affect the peak positions of Cl(Q,lo) (see, 
for instance, Fig. (|27|l for the aluminum case). 

Inelastic X-ray Scattering, with lower energy resolu- 
tion (a few hundreds of meV) and in a broader energy 
range (up to a few eV), allows to study electronic exci- 
tations (plasmons). In this case the plasma oscillation 
is brought about by free electrons, while the background 
is c onstituted by the ionic n etwork. Some recent stud- 
ies l|Burns et flit 119991 12001 point out a decrease of the 
plasmon dispersion at low metal concentrations, which, 
in turn, is ascribed to the failure of the RPA approxima- 
tion. 



V. SUMMARY AND PERSPECTIVES 

In this section we will try to summarize the scenario 
arising from the measurements reported so far. In respect 
to the collective properties, it seems useful to discuss the 
different dynamical regimes probed at different wavevec- 
tors and frequencies. Although the two domains are in 
principle independent (as testified, for instance, by the 
two separate generalization of the classical hydrodynam- 
ics introduced in sec. III.F|I . the existence of a disper- 
sion relation ultimately allows one to think in terms of 
wavelenghts only. In our point of view, though precise 
boundaries can not be traced, one can identify in liquid 
metals the following dynamical regimes: 



1 -Hydrodynamics 



• The hydrodynamic, Q — > limit, that, in liq- 
uid metals, basically means Q < 0.1 nm . In 
this region simple hydrodynamic treatment based 
on three microscopic dynamical variables (density, 
current, energy) provides an exhaustive description 
of the main features. Although not accessible by 
neutron and X-ray spectroscopic techniques, the 
hydrodynamic predictions are in very good agree- 
ment with acoustic measurements of sound velocity 
and attenuation properties. Moreover, the strict 
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FIG. 38 Sketch of the different dynamical regimes on decreas- 
ing the wavelength. 



analogy of long wavelength fluctuations in con- 
ductive and ordinary liquids (accessible via light 
scattering) corroborates the validity of the sim- 
ple hydrodynamic theory. This region is charac- 
terized by adiabatic sound propagation, and the 
whole dynamical features are ruled by macroscopic 
transport parameters (viscosity, thermal diffusivity, 
specific heats). In this regime, other approaches 
tailored for conductive fluids such as plasma os- 
cillation theories provides alternative descriptions, 
which turns out to be increasingly accurate for low 
r s systems. 

An (hypothetic) isothermal region, which should 
be observable in the 0.2 < Q < 3 nm" 1 momen- 
tum range. Upon increasing the wavevector out- 
side the strict hydrodynamic limit, indeed, the life- 
time of the entropy fluctuations becomes increas- 
ingly shorter (r = (jDtQ 2 )~ 1 )- Since the fre- 
quency of the corresponding density fluctuation in- 
creases almost linearly {u w cQ), one should expect 
a transition at Q* w c/Dt- Given the typical ther- 
mal diffusivity values of liquid metals this crossover 
should be located at a few fractions of nm -1 . Since 
with both INS and IXS one normally access mo- 
mentum transfers above 1 nm -1 , this region can 
be hardly explored, and no direct and convincing 
indication for its existence are avai lable at prese nt. 
Some old INS data on liquid lead l)Faberl . ll972f) as 
well as a more recent experiment on liquid nickel 
(B ermeio et all l2000j) . however, seems to suggest 
indication for such existence. The case of Nickel, 
in particular, is an interesting one and would de- 
serve new investigations since in this system the 
specific heat ratio is particularly high and there- 
fore the difference in sound velocities between an 
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ordinary adiabatic regime and an isothermal one 
would be of the order of 40%. It is worth to stress 
that the prediction for the existence of an isother- 
mal regime poses on: i) a negligible Q dependence 
of the thermodynamic quantities and ii) the va- 
lidity of a one component effective description in 
which the effective thermal diffusivity is well de- 
scribed by the sum of the electronic and ionic con- 
tributions. Both this points have been recently 
discussed in the analysis of alkali metals and liq- 
uid a luminium IXS spectra llSco pigno a nd Ruoccol 
l2004USingh and Tankeshwaii 120031 l2004h . 

• A generalized hydrodynamic region, probed above 
Q « 3 nm , in which the frequency- wavevector 
dependence of the transport properties heavily af- 
fects the sound mode. The upper limit of valid- 
ity of this region is rather system dependent, nor- 
mally higher for alkali metals (say up to 0.7Q m )- 
The natural framework is here the memory func- 
tion formalism, which can be developed at differ- 
ent levels of accu racy, ranging fr om the celebrated 
Lovesey's model ijLovesevl Il97l|) (accounting for a 
single average relaxation time for the second or- 
der memory function) up to more refined mem- 
ory function models based on multiple relaxation 
phenomena whi ch, firstly introduced for Lennard- 
Jones systems ijLevesaue et all Il973jh have been 
more recently adapted and t ested against IXS in- 
vestigations of liquid metals ( [ Monaco et all . l2004t 
IScopjgno er^] . l2000bicl l2002b) and numerical sim- 
ulations of model und ercooled and glassy alkali 
l)Scopigno et all\2002^ . 

In this respect, it is worth to point out how the high 
points of a memory function approach are not solely 
in a better agreement with experimental data, 
which, in general, heavily depends on the quality 
of the experimental data and, of course, on the 
numbe r of the model parameters ijScopigno et all 
I2000cl) . In most cases, indeed, simplified phe- 
nomenological models such as the damped har- 
monic oscillator and a lorenzian function for the 
inelastic and quasielastic feature, respectively, pro- 
vide satisfactory agreements. 

On the contrary, the memory function framework 
allows to grasp an insight behind the mechanisms 
ruling the relaxation dynamics, extracting relevant 
information about quantities which are not directly 
related to any spectral features, such as the Q de- 
pendencies of the relaxation time(s) and of the lon- 
gitudinal viscosity. 

Following the prescriptions illustrated in section 
III.F.2I indeed, the longitudinal viscosity is related 
to the total area under the viscous contribution to 
the memory function. In Fig. (|39|l . for instance, we 
report the (generalized) longitudinal viscosities ex- 
tracted in this way for several investigated systems. 
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FIG. 39 Generalized longitudinal viscosities for a collection 
of liquid metals, extracted by the experimental data through- 
out a memory function approach. The long wavelenght limit, 
determined by ultrasonic measurements, is also reported for 
K, Na, Ga (bottom to top). 



As can be seen, the low Q extrapolation of the ex- 
perimental values compares well with independent 
acoustic determinations (when available), but also 
allows to determined the Q generalization of this 
important transport properties, which can be di- 
rectly estimated only by numerical simulations ap- 
proaches. 

In this Q region, hard sphere-based theories pro- 
vide alternative descri ptions in terms of ex t ended 
hydrodynamic models ijKamgar-Parsi et~all Il987|) 
but they still miss a convincing explanation of one 
of the most important points: the physical interpre- 
tation beyond the relaxation of the sound mode, 
which is now a firmly established evidence sup- 
ported by uncountable experimental investigations. 

This latter aspect, which usually emerges in terms 
of a speed up of the sound velocity taking place in 
the 1-9-10 nm" 1 region, is one of the most inter- 
esting aspects which is still lacking an explanation. 
Broadly speaking, it is always referred to as a shear 
relaxation, but the ultimate nature of the involved 
physical processes still have to be clarified. 

Actually, Mod e Coupling Theory 

(jErnst and Dorfmarl Il975(l provides a description 
of the acoustic dispersion curve in te rms o f even 
powers of Q (|Ernst and Dorfmanl Il975j) . but 
its interpretation is restricted to monoatomic 
systems in the liquid phase, while the observed 
phenomenology seems to be a more general feature 
of the disordered systems. 

Interestingly, indeed, a qualitatively similar phe- 
nomenon (see Fig. I40|) has also been re- 
ported numerically in fused silica l)Horbach et all 
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FIG. 40 Soun d dispersion in severa l kind of liquids. Al- 
kali m etals (Li flScopigno et all\2QCtCl(t ) and K ijMonaco et all 
l2004fl . dashed line is the iso thermal sound veloci ty); Hydro- 
gen bonding s ystems (Water j Monac o et aZ.lll999t) and hydro- 
gen Fluoride (Aneelini et al, 2002)) low and high frequency 
sound velocities are also ind icated by the dotted lines); liquid 
neon llCunsolo et aZ.ll200lD : adiabatic (circles and line) ap- 
parent (open triangles) and high frequency (open diamonds) 
sound velocities ; liqui d silica, molecular dynamics simulations 
dHorbach et ffll.Ll200lD . in this case also the transverse branch 
is reported. 



l200ll)) and experimentally in Lennard jones flu- 
ids llCunsolo et HI l200ll) and h yd rogen bond- 
ing s ystems ( ijAngelini et aH l2002t iMonaco et all 
1999)). While in this latter case the positive dis- 
persion has been shown to be an activated pro- 
cess, related to the structural relaxation, a differ- 
ent scenario seems to characterize the other sys- 
tems. Interestingly, indeed, (see Fig. I41|l the 
same behavior of the sound velocity also appear in 
glasses either in experiment al IXS measur e ments 
in g-Se and vitreous silica ijRuzicka et all l2004t 
IScopigno efall . l2004|) . in numerical works on model 
glass es of Lennard J ones systems (iRuocco et all 
12000^ . vitre ous silica llHorbach et ^J . l200l[T and al- 
kali metals l|Scopigno et all l2002c|) and, finally, in 
theoretical calculations for an hard sphere glass 
i|Gotze and Mavrl fcOQO). In the case of alkali met- 
als, in particular, it has been shown how the speed 
up of the sound velocity persists upon cooling well 
below the glass transition, thus ruling out the pos- 
sible role of the structural relaxation in this effect. 
The presence of positive dispersion at THz frequen- 
cies in glasses, quantitatively comparable to the one 
observed in liquids, therefore, challenges the inter- 
pretation of the positive dispersion in terms of a 
transition from a liquid like to a solid like behav- 
ior, an effect which seems to be quantitatively neg- 
ligible (with the remarkable exception of hydrogen 
bonding systems). Accepting a description of the 
collective dynamics proceedings over two distinct 
viscous relaxations, therefore, the ultimate respon- 
sible for the bending up of the dispersion curve 



seems to be the faster of the two observed processes, 
which turns out to have a mild temperature depen- 
dence. In this respect, the physical nature of this 
faster process calls for deeper understanding. It is 
worth to point out however, that in Lennard Jones 
systems the positive dispersion is recovered within 
an harmonic description of the dynamic structure 
factor in terms of eigenvalues and eigenvectors, a 
result which relate the positive dispersion to the 
properties of the dynamical matrix and, ultimately, 
to the topological disorder of the inherent equilib- 
rium position, being the interaction potentials com- 
parable in glasses and crystals. 




FIG. 41 Sound velocities in several glass y systems. Hard 
spher es glass, exact solution within MCT llGotze and Mavil 
2000); Metallic glass obtained quenching a model system 
intera cting via Price- Tosi pseudopotential llScopigno et all 
l2002cT ). the crystalline counterpart is a lso reported; Lennar d 
Jones glass obtained in a similar way (iR.uocco et all l2000t) . 
reported with the low and high frequency s ound velocities; 
Si02, MD as in Fie l4()l but in the glassy state ijHorbach et all 
l200lf) a nd experimentally de termined by means of IXS scat- 
tering ijRuzicka et all I2004T) . In both case s the transverse 
branc h is also reported; glassy selenium llScopigno et all 
I2004T) . again the sound velocity exceeds the adiabatic value. 

• A kinetic regime, valid from around Q m , up to 
a few oscillations of the structure factor. Here 
the hard sphere description provides remarkably 
accurate predictions, in terms of a quasielastic 
"extended heat mode" whose linewidth is de- 
scribed in terms of the Enskog's diffusion coef- 
ficient and of an equivalent hard sphere diam- 
eter. The extent of validity of kinetic theory 
in this momentum range has been widely tested 
(ICohen et a/.Ul984ll987ljde Schepper et a/lll984t 
iKamgar-Parsi et all Il987(l in several hard sphere- 
like systems (like alkali metals and lennard-jones 
fluids), and it also apply in colloidal systems. It 
would be interesting to challenge such theory in 
less simple liquid metals. Despite the success of 
hard sphere model, there is still an obscure point 
concerning the real origin of such extended heat 
mode: while in the low Q limit it tends to the en- 
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tropy fluctuation mode, indeed, at finite Q's it cer- 
tainly involves mass diffusion processes. Looking 
at things from the constant energy point of view, 
umklapp modes resembling crystalline phonons in 
Brillouin zones high er than the first seems to be still 
poorly under s tood llCock ing and EgelstaflL [l965a; 

[ Dorner et all 19651: iRandolph and Sinew! 11966: 
Scopieno et adl2001b|) ~ 



• An high Q region, probed as soon as the Van Hove 
distinct function vanishes. In this limit, both IXS 
and INS experiments provide the same information, 
about the atomic motion over timescales shorter 
than the interparticle collision time. An almost ex- 
act description for this regime is available, which 
can also account for quantum aspects such as re- 
coil energy an quantum corrections to the spectral 
moments. 
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The single particle dynamics seems to be better under- 
stood if compared to the collective motion, at least by a 
coarse grained point of view of an hydrodynamic diffu- 
sive mode with finite Q corrections, evolving towards a 
ballistic regime. This is probably due to the intrinsic 
difficulty of isolating a (wide enough) constant Q coher- 
ent energy spectrum from an INS experiment. Since the 
study of the strictly coherent spectrum became possible 
only in the last decade, indeed, a lot of efforts have been 
devoted in the past to the single particle case. Though 
all the approaches described in this review describe, on 
average, equally well the available experimental data, a 
memory function formalism paralleling the one for the 
collective case could provide the route to relate the sin- 
gle particle motion and the collective dynamics in the 
microscopic regime. 

Raising the level of detail of the description of the 
single particle dynamics, however, a major experimen- 
tal challenge seems to be the identification of the dif- 
ferent processes giving rise to the quasielastic incoher- 
ent scattering. Recent INS results, indeed, suggest the 
presence of two distinct physical mechanisms, active 
over different timescales , underlying the diffusive mo- 
tion ijBadval et q?l 120031: iBove et all 12001 l2002al) . The 
combined presence of coherent and incoherent scattering, 
however, makes such identification still unclear although, 
in principle, the IXS signal might be used to subtract 
the coherent contribution from the INS spectra, thus ex- 
tracting the purely incoherent dynamics. In this respect, 
the synergy of combined IXS and INS studies on a same 
sample seems us imperative and could help to shed light 
on this point. The IXS signal indeed, might be used to 
subtract the coherent contribution from the INS spectra, 
thus extracting the purely incoherent dynamics. 
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